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We perform a general study of the thermodynamic properties of static electrically charged black 
hole solutions of nonlinear electrodynamics minimally coupled to gravitation in three space dimen- 
sions. The Lagrangian densities governing the dynamics of these models in fiat space are defined 
as arbitrary functions of the gauge field invariants, constrained by some requirements for physical 
admissibility. The exhaustive classification of these theories in flat space, in terms of the behaviour 
of the Lagrangian densities in vacuum and on the boundary of their domain of definition, defines 
twelve families of admissible models and allows for a complete characterization of the associated 
sets of gravitating electrostatic spherically symmetric solutions by their central and asymptotic be- 
haviours. We focus on nine of these families, which support asymptotically Schwarzschild-like black 
hole configurations, for which the thermodynamic analysis is possible and pertinent. In this way, the 
thermodynamic laws are extended to these cases and the generic behaviours of the relevant state 
variables for the associated sets of black hole solutions are classified and thoroughly analyzed in 
terms of the aforementioned boundary properties of the Lagrangians. Moreover, we find universal 
scaling laws (which hold for any admissible model) running the thermodynamic variables with the 
electric charge and the horizon radius. The scale transformations form a one-parameter multiplica- 
tive group, leading to universal "renormalization group"-like first-order differential equations. The 
beams of characteristics of these equations generate the full set of black hole states associated to 
any of these gravitating nonlinear electrodynamics. Moreover the application of the scaling laws 
allows to find a universal finite relation between the thermodynamic variables, which is seen as a 
generalized Smarr law. Some particular well known (and also other new) models are analyzed as 
illustrative examples of these procedures. 

PACS numbers: 04.70. Bw, 04.70.Dy, ll.10.Lm 



I. INTRODUCTION 

In General Relativity the theorems on singularities, on 
the uniqueness of the black hole solutions of the Ein- 
stein vacuum field equations, and the Cosmic Censor- 
ship conjecture [l[ lead to the strong belief that the Kerr 
spacetime describes accurately the outcome of the grav- 
itational collapse of rotating electrically neutral matter 
distributions. When such distributions support a non- 
vanishing total electric charge, similar arguments, now 
involving the minimally coupled Einstein-Maxwell the- 
ory, lead to Kerr-Newman black hole configurations as 
the final state of the collapse 2] . 

However, there are several arguments justifying the 
study of other BH configurations @, which are solu- 
tions of more complex gravitating NEDs, as possible final 
states of the gravitational collapse of electrically charged 
matter. A first argument comes from the low energy 
regime of string theory. If one accepts the fundamental 
character of this theory, then the dynamics of the emerg- 
ing gauge fields in this regime is governed by nonminimal 
effective Lagrangians |4J. In this context, some nonlinear 
electrodynamics, as the Born-Infeld model @, are to be 
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seen as more "fundamental" in describing the dynamics 
of the electromagnetic fields than the usual Maxwell the- 
ory, which is regarded, within this picture, as a weak field 
approximation. 

A second argument concerns the quantum vacuum ef- 
fects on the electromagnetic field, which are neglected 
in the classical Einstein-Maxwell theory. Now accepting 
the fundamental character of the Maxwell theory and its 
quantum formulation, the corrections of the Dirac vac- 
uum on the dynamics of electromagnetic fields in flat 
space can be described, at a phenomenological classical 
level, through nonlinear effective Lagrangians that de- 
fine NEDs, the first historical example being the EH La- 
grangian Q. It accounts for the nonlinear effects of the 
Dirac vacuum on the low energy electromagnetic wave 
propagation, calculated to lowest order in the pertur- 
bative expansion. When higher order operators are in- 
cluded in the perturbative calculation we are led to a 
sequence of effective Lagrangians which take the form of 
polynomials in the field invariants, arranged as an ex- 
pansion in operators of increasing dimensions Q. These 
effective Lagrangians exhibit finite-energy elementary so- 
lutions, even though the bare (Maxwell) Lagrangian does 
not [8j. When minimally coupled to gravity such La- 
grangians describe satisfactorily the vacuum effects in 
the slowly varying curvature regions [9]. Disregarding 
the rotational degrees of freedom, which are not essential 
for the present considerations, the elementary solutions 
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of any physically admissible NED (a concept to be de- 
fined below) minimally coupled to gravity, are asymptot- 
ically flat, exhibit a curvature singularity at the center 
and their electrostatic field strengths are monotonically 
decreasing functions of the radius [Tol . [TTI j . Near the sin- 
gularity this effective approach fails, and the description 
of the vacuum effects would require a consistent quantum 
treatment of both the electromagnetic and gravitational 
fields, which is not available at present. But in moving 
away from the singularity the curvature stabilizes very 
quickly and becomes slowly varying in the characteris- 
tic distances of the vacuum screening. Consequently, the 
effective approach remains accurate everywhere in the 
manifold, excepted in a small region close to the singular- 
ity. Therefore in describing the collapsed state of nonro- 
tating electrically charged matter, the BH configurations 
associated to some NEDs could be more "realistic" than 
the Reissner-Nordstrom one. 

A general analysis of the gravitational configura- 
tions and geometric structures of the elementary so- 
lutions associated to the whole class of admissible 
NEDs, minimally coupled to gravity (including some well 
known models such as Maxwell, Born-Infeld or Euler- 
Heisenberg, as particular representative examples) has 
been performed in Refs.[lCj, Nevertheless, the issue 
of the thermodynamic properties and laws of the corre- 
sponding BH solutions was not considered in those works. 
Such a study is the main object of this paper. 

The thermodynamic analogies for asymptotically flat 
BH solutions of the Einstein equations were established 
in the early seventies for the cases of the Kerr-Newman 
BH and its nonrotating and uncharged limits 12]. The 
subsequent discovery of Hawking's radiation [1 31 ] turned 
these analogies into actual laws of the thermodynamics of 
these configurations. Since then, an important amount of 
literature has dealt with this issue in depth. In this way 
many papers have analyzed the thermodynamics of the 
Reissner-Nordstrom black hole (see e.g. [3] and refer- 
ences therein, also [HI). Let us mention some extensions 
and developments going beyond the analysis of this solu- 
tion: i) the validity of the thermodynamic laws has been 
established for asymptotically flat BH solutions of sev- 
eral classes of field theories minimally coupled to gravity, 
including NEDs supporting asymptotically coulombian 
elementary solutions (l6| . In this way the thermody- 
namic analysis has been performed for asymptotically 
flat BH solutions associated to some particular NEDs, 
as the gravitating BI model [l7| or, more recently, the 
family of NED models defined through Lagrangian den- 
sities which are constructed as powers of the Maxwell one 
[l8j |: ii) this analysis has been also performed for BH so- 
lutions of the gravitating Born-Infeld NED in asymptoti- 
cally cosmological backgrounds, as (anti-)de Sitter spaces 
in (3+1) dimensions |l9l [20|, in higher dimensions [2ll[22| 
and in (2 + 1) dimensions [23J, extending the Hawking- 
Page results of the uncharged case [24| and motivated 
by the AdS/CFT correspondence [25| : iii) several exten- 
sions of these solutions, resulting from the minimal cou- 



pling of NEDs to gravitational dynamics defined beyond 
the Einstein-Hilbert action, have been considered, as in 
the cases of Gauss-Bonnet and Lovelock (2(1, f(R) H3 



or scalar-tensor theories [28j ; iv) statistical computations 
of the entropy for the extreme and near extreme five di- 
mensional RN black holes have been performed in the 
framework of string theory [29j . 

However, most of the aforementioned works deal with 
some particular models, while a comprehensive study of 
the thermodynamic properties of BH configurations as- 
sociated to general gravitating NEDs is possible, but 
still lacking. The aim of this paper is to provide such 
a study, by analyzing the full set of NEDs minimally 
coupled to gravitation in (3 + 1) dimensions and defined 
by Lagrangian densities which are arbitrary functions of 
the two field invariants, but restricted by some consis- 
tency requirements (such as the positivity of the energy, 
the regularity of the Lagrangian and the parity invari- 
ance) leading to what we shall call admissible theories. 
In Ref.[ll[ the gravitating admissible NEDs were classi- 
fied into two sets, according to the form of the asymp- 
totic behaviours of their elementary solutions, which are 
always asymptotically flat, but can approach flatness as 
the Schwarzschild field (normal case) or slower than it 
(anomalous case). In the anomalous case some thermo- 
dynamic variables (such as the ADM mass) cannot be de- 
fined and the thermodynamic analysis is not possible, at 
least in the usual way. Consequently, in the present work 
we shall deal only with the analysis of the thermodynam- 
ics of asymptotically normal BH solutions supported by 
the corresponding families of G-NEDs. 

In section [H] we summarize results of Ref . [8| regarding 
the elementary solutions of NEDs in flat space, which 
are necessary for the present study. In particular, we 
introduce the classification of these theories that, when 
generalized to the case of G-NEDs, defines the different 
classes of possible BH solutions. 

In section lnTl we briefly describe the geometrical struc- 
tures of the BHs associated to each family, which are 
necessary for their thermodynamical analysis. A detailed 
report of the contents of this section has been published 
in Ref. PH. 

Section IIVI is devoted to the thermodynamical analy- 
sis itself. First, we prove the validity of the expression 
of the first law for any admissible, asymptotically nor- 
mal, gravitating NED, in particular for those with non 
Maxwellian weak field limit, for which the theorems of 
Ref.[l6| do not immediately apply. Next we analyze the 
behaviour of the main thermodynamical variables and 
specific heats for the BH configurations associated to the 
different families. In particular the study of the thermal 
properties of the extreme BH configurations suggests the 
possible presence of van der Waals-like first-order phase 
transitions between these states for several families. 

In section [V] we find some universal relations between 
the thermodynamic variables of the asymptotically- 
normal BH states associted to admissible G-NEDs. First, 
we obtain a universal expression for the generalized 
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Smarr law, valid for all such BH solutions. This law 
reduces to the usual Smarr formula in the RN case and 
to the generalized expressions obtained for the few par- 
ticular examples of G-NEDs studied in the literature. 
Next, we state universal scaling laws satisfied by the 
functions relating the thermodynamic variables for any 
G-NED. We make a preliminary discussion about some 
implications of these laws. In particular, we show the 
existence of a one-parameter multiplicative group struc- 
ture underlying the scaling transformations and obtain 
the associated first-order differential equations ("renor- 
malization group" -like equations) relating the derivatives 
of the state variables with respect to the charge and the 
horizon radius. The characteristic curves of these equa- 
tions define the running of the thermodynamic variables 
and generate the "equations of state" for the sets of BH 
solutions associated to any gravitating NED. 

Finally, we conclude in section lVTl with some comments 
and perspectives. 



dX ~ dY 2 ' 



(3) 



(for X — E 2 (r > 0),Y — 0) which is necessary and 
sufficient for the linear stability of the ESS solutions, as 
shown in Ref.Q. 

The Euler field equations associated to the La- 
grangians ([!]) read 



dX dY 



= 0, 



(4) 



and the corresponding energy-momentum tensors take 
the form 



dip 



dip 



8X 



dY 



(•5) 



II. GENERAL RESULTS ON NEDS IN FLAT 
SPACE-TIME 



In looking for ESS solutions of the field equations (E(r) = 
E(r)f/r, H = 0) we are led to the first integral 



Let us summarize the results concerning the complete 
characterization and classification of all the families of 
admissible NEDs through the properties of their ESS so- 
lutions [8] , which will be necessary for the analysis of the 
gravitating problem. The dynamics of these fields in flat 
four-dimensional space-time are governed by Lagrangian 
densities which are single-branched functions 



ip(X,Y), 



(1) 



of the two gauge field invariants X = — \F llill F tiU = 

E 2 - H 2 and Y = -\F^F*^ V = 2E ■ H, where 
Ffiu = dfj,A v — dyA^ is the field strength tensor of the 
vector potential A M , F*^ v = \^ va ^F a& its dual and 
E and H the electric and magnetic fields, respectively. 
For physical consistence of the corresponding theories 
these functions must be defined in a open and con- 
nected domain of the X — Y plane, including the vacuum, 
and be restricted by the conditions of parity invariance 
(ip{X, Y) = ip(X, —Y)) and positivity of the energy. The 
latter condition reads M 



> VX 2 + Y 2 + X 



dip 
dX 



Y 



dip 
dY 



cp(X,Y)>0, (2) 



with vanishing vacuum energy (ip(0, 0) = 0). Moreover, 
the ESS solutions E(r) must be single-branched, defined 
and regular for r > 0. This requires ip(X,0) to be a 
strictly monotonically increasing function of X in the 
range of values X = E 2 (r > 0) {dip/dX \y=o is required 
to be continuous and strictly positive there). The above 
set of requirements define what we have called admissi- 
ble NEDs, the only ones that we shall consider here. In 
addition one can assume ip(X, Y) to satisfy the condition 



i dip 
dX 



Y=0 



E(r) = Q, 



(6) 



where the integration constant Q is identified with the 
point-like charge source of the solution. Once the form of 
ip(X,Y) is specified Eq.© gives in implicit form the ex- 
pression of the ESS field E(r, Q) as a function of rj \ 
in such a way that it scales as 



\E(r,Q)\ = \E(R,Q = l)\, 



(7) 



(R 



As we shall see below, this scaling law in- 



duces similar laws in the different thermodynamic func- 
tions describing the systems of BH solutions of any gravi- 
tating NED. We shall call "characteristic" the Q = 1 con- 
figurations, from which the whole thermodynamic struc- 
ture of the set of BH solutions associated to a given G- 
NED is determined through these scaling transformations 
(see sections [TV] and [V] below) . 

The positivity of the derivative of the Lagrangian func- 
tion in Eq.© allows to restrict the analysis to the case 
(E > 0, Q > 0) without loss of generality. Assuming 
behaviours for the ESS solutions around the center and 
asymptotically of the form 



E(r -> 0,Q) - vi{Q)r p ; E{r -> oo, Q) - ^{Q)r q , 

the admissibility conditions restrict the values of the ex- 
ponents to p < and q < 0. At the center the fields 
diverge for p < and behave as 



E(r -> 0) ~ a - b(Q)r° 



(9) 
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for p = 0. In the latter case the parameter a (the maxi- 
mum field strength) and the exponent a > are universal 
constants for a given model, whereas the coefficient b(Q) 
is related to the charge of each particular solution as 



KQ)Q a,2 = Hm 

X— >a- 



(a - VI) 



er/2 



(10) 



&o = b(Q = 1) being also a universal constant of the 
model. Asymptotically E(r —> oo) vanishes for the ESS 
solutions of any admissible model. 

The corresponding behaviours of the Lagrangian den- 
sities around these values of the fields (X = E 2 (r)) are 
given by the expression 



tp(X, Y = 0) ~ ctiX^ 



(11) 



where oii,{i = 1,2) are positive constants which are re- 
lated with the coefficients in Eq.®. For X — > oo (and 
p < 0) we have 



for the energy density distributions of the ESS solu- 
tions. As easily seen, the spatial integrals of these func- 
tions around the center of the ESS solutions converge if 
— 1 < P < and diverge if p < — 1 . This allows a classifi- 
cation of the admissible NEDs in terms of the central be- 
haviour of their ESS solutions, leading to three families. 
We have called Al the first family, corresponding to the 
range of values — 1 < p < 0, for which the ESS solutions 
diverge as r — > 0, but the integrals of energy converge 
there. For these models the Lagrangian functions be- 
have, when X — > oo, as in Eq. fTTj) with 7 > 3/2. We have 
called A2 the second family, for which p = 0, correspond- 
ing to the NEDs supporting ESS solutions which are fi- 
nite at the center. Obviously, the integrals of energy con- 
verge there for these solutions (see Eqs. ([T4^) - (IT51) '). The 
third family corresponds to p < —1 and supports ESS 
solutions diverging at the center and whose integrals of 
energy diverge also there (we have called these models 
UVD, acronym for ultraviolet divergent). When the in- 
tegral of energy is convergent at the center (cases Al and 
A2) the "internal energy function" 



MQ) 



7if*i 

Q 



p/2 



; 71 



whereas for X — > we have 



q/2 
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p-2 
2p ' 



g-2 
2q ■ 



(12) 



(13) 



In both cases 7, > 1/2. If p = the Lagrangian densities 
behave around the center (X = E 2 (r = 0) = a 2 , Y = 0) 



[p (X,Y = 0)~ 2 -^(a-Vx)^+A, 
if (7 7^ 2. If a = 2 this behaviour becomes 



<p(X, Y = 0) ~ -6 Ha ~ VI) + A 



(14) 



(15) 



In these equations A = ip(a 2 ,Y = 0) < 00 if a > 2, 
and it is a universal constant of the model. If a < 2 the 
value of A can be calculated from the explicit expression 
of the Lagrangian density (see Eqs. (l8Tj) and (1821) below). 
Consequently, for a < 2 the Lagrangians exhibit a ver- 



tical asymptote at X 



For a > 2 they take finite 



values there, with divergent slope (the Born-Infeld model 
5] belongs to this case with a = 4). 
Equation (J5J) leads to the expression 



- T° - 2^ 



Y=0 



E 2 (r,Q)-cp(X,Y) 



2QE - r 2 y 
(16) 



e m (r,Q)=4ir / R 2 p(R,Q)dR, 



(17) 



giving the energy of the ESS field of a point-like charge 
Q inside a sphere of radius r, is well defined. 

The asymptotic behaviour of the ESS solutions al- 
lows for a second classification of the admissible NEDs. 
For large r the integrals of energy converge asymptot- 
ically if q < —1 and diverge if —1 < q < 0. Thus 
we have defined four classes of models differing by this 
behaviour: Bl models (when —2 < q < —1) whose 
ESS solutions are asymptotically damped slower than 
the Coulomb field and the integrals of energy converge 
asymptotically; B2 models (when q — —2) whose ESS 
solutions behave asymptotically as the Coulomb field; 
B3 models (when q < —2) whose ESS solutions vanish 
asymptotically faster than the Coulomb field; IRD mod- 
els (acronym for infrared divergent) when — 1 < q < 0, 
whose ESS solutions are asymptotically damped slower 
than the Coulomb field and their integrals of energy di- 
verge asymptotically. When the integral of energy con- 
verges asymptotically (cases Bl, B2 and B3) the "ex- 
ternal energy function" 



/>oo 

: (r,Q)=47r/ R 2 p(R,Q)dR, 

J r 



(18) 



giving the energy of the field outside the sphere of radius 
r, is well defined. 

The combination of these central 
haviours allows for the exhaustive 
admissible NEDs in twelve classes, 
nations of Al and A2 with Bl, B2 
support finite-energy ESS solutions, 
terpreted as nontopological solitons 



and asymptotic be- 
classification of the 
Six of them (combi- 
and B3 behaviours) 

which are to be in- 
if the necessary and 
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sufficient stability condition <|3j> is fulfilled by the corre- 
sponding Lagrangian densities in the domain of definition 
of these solutions Q . For these six families the conver- 
gence of the integral of energy in the whole space leads 
to the relations 



e(Q) = e m (oo, Q) = e ex {0, Q) = s in (r, Q) + s ex (r, Q) 

(19) 

for the energy functions associated to the ESS solutions 
of charge Q. Using Eqs.dH]) and (fTl))) it can be easily seen 
that these energies, when they are well defined, scale as 



IP (X, Y = 0) 




e(Q) = Q^ 2 e{Q = 1); 



e ln (r,Q) = Q 3/2 e l 
£ e *(r,Q) = Q 3/2 £ e 



,Q = l 



Q = \ 



where the total energy of the characteristic configuration 
e(Q = 1) is a universal constant for a given model. By 
using the last expression of g in Eq. (|16[) for the calcu- 
lation of (fT8|) as r —> 0, after some partial integrations 
and taking account the first integral (j6|), as well as the 
central and asymptotic behaviours of the fields in these 
finite-energy cases, we are led to the useful expression 



Figure 1. Qualitative behaviour of the admissible Lagrangian 
densities tp(X, Y = 0) i) for small ESS fields {E 2 = X ~ 0, 
corresponding to the B and IRD asymptotic behaviours of 
( 20 ) the ESS solutions), ii) for large ESS fields {E 2 = X -> oo, cor- 
responding to the Al and UVD central-field behaviours) and 
iii) for finite maximum field-strength models (X = = 
a 2 ), corresponding to the A2 central-field behaviour. The La- 
grangian densities behave as X 72 around X = (in the B and 
IRD cases) and as X 71 around X — > oo (in the Al and UVD 
cases) where the 7; constants are related to the central and 
asymptotic behaviours of the ESS solutions through Eqs. (|12[) 
and (|13[) respectively. In the A2 cases the Lagrangian den- 



e{Q) = 



16ttQ 



E(R,Q)dR, 



(21) 



sity exhibits a vertical asymptote at X 



(if a < 2) or 



takes a finite value with divergent slope there (if a > 2). In 
the intermediate range of X values, matching the central and 
asymptotic regions, tp(X, Y = 0) must be strictly monotoni- 
cally increasing, for admissibility. 



relating the total energy of the soliton with the integral 
of the field. 

The ESS solutions of the remaining six cases (combi- 
nations including UVD, IRD or both behaviours) are 
energy divergent. When coupled to gravity all the models 
of the twelve families lead to consistent gravitating field 
theories. However, we shall see that the families involv- 
ing the IRD behaviour fail to allow for a thermodynamic 
analysis in the usual way. 

In figure Q] we have plotted the profiles of the La- 
grangian density functions of the admissible NEDs (for 
Y = 0) around the vacuum and on the boundary of their 
domain of definition. 



As usually, the action corresponding to the minimal 
coupling of gravity and any NED reads 



S — Sc4 + SnED 



R 



16nG 



(22) 



g and R being the metric determinant and the curvature 
scalar, respectively. The associated Einstein equations 
for the G-ESS solutions can be written in an adapted 
Schwarzschild-like coordinate system, with an interval of 
the form 



III. GENERAL RESULTS ON GRAVITATING 
NEDS 

The full set of admissible NEDs minimally coupled to 
gravity, and the structure of their G-ESS solutions have 
been classified and extensively analyzed in Refs.jlOl. Illj. 
However, as mentioned in section HI this analysis con- 
cerned mainly the geometric structure of these solutions 
and excluded the study of their thermodynamic proper- 
ties. Before tackling this thermodynamic problem let us 
summarize some pertinent results of these references. 



ds 2 =g(r)dt 2 -^-- r 2 dtf, (23) 
9{r) 

where d£! 2 = d9 2 + sin 2 '(9)dip 2 is the angular metric and 
5oo = g{f) is the only metric function to be determined. 
The simplicity of the form of this interval is a consequence 
of the Einstein equations together with the equality of 
the mixed components of the energy-momentum tensor, 
Xq = T/, which remains the same as in flat space in 
spherical coordinates. These equations reduce to 
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d 
dr 



{rg(r) - r) = -8^r 2 T ° = -87rr 2 ( 2^E 2 s 



Q c = 



Wira 



(29) 



dr 2 



(rg(r)) — — WirrT 2 — 16nrtp, 



(24) 



whose compatibility can be easily established. The gen- 
eral solution of this system reads 



g(r, Q, C) = 1 + — - — / r 2 T °(r, Q)dr, 
r r J 



(25) 



where C is an arbitrary integration constant. The char- 
acterization and classification of the admissible NEDs 
in flat space, summarized in section [Til allows for the 
complete characterization and classification of the cor- 
responding G-NEDs. The asymptotic behaviour of the 
metric function is dominated by the behaviour at large 
r of the ESS solutions in flat space and is independent 
of their properties at finite r. Indeed, when the asymp- 
totic behaviour of the ESS solutions in flat space belongs 
to the cases Bl, B2 or B3 the metric function can be 
written as 



g(r, Q,M) = 1 1 , (26) 

r r 

where e ex (r, Q) is the external energy function, defined 
by (p~5|) . and the constant M is identified as the ADM 
mass. Owing to the asymptotic behaviour of the external 
energy function in these cases (~ r q+1 ,q < — 1), the last 
term in ((26)) can be neglected at large r and the metric 
function behaves asymptotically as the Schwarzschild one 
[l0| . The combinations of these asymptotic B-cases with 
the central-field Al and A2 cases, which support finite- 
energy ESS solutions in fiat space, admit an alternative 
expression for the metric function in terms of the internal 
energy function, which reads 



g(r,Q,M) = l , (27) 

r r 

where the constant U is now related to the gravitational 
mass and the electrostatic soliton energy in flat space 
through 



U = M-s(Q). 



(28) 



This constant may be roughly interpreted as the gravita- 
tional binding energy of the configuration plus the proper 
energy of the point-like source (see, however, the discus- 
sion about this kind of interpretations and the associ- 
ated difficulties in Rcf.(3fJ and the pertinent references 
therein). In reference [l(| we have called "critical" the 
configurations with U = 0. As we shall see in the next 
section in the A2 cases with the particular value of the 
charge 



these configurations correspond to extreme black points 
[31], which exhibit a rich and peculiar thermodynami- 
cal behaviour. With Q < Q c they correspond to non- 
extreme BPs. In A2 cases with Q > Q c (or for Al cases 
with any Q), they correspond to BH configurations with 
an external event horizon and a vanishing-radius Cauchy 
inner horizon. 

For models belonging asymptotically to B-cases and 
with UVD central-field behaviour (e.g. Maxwell La- 
grangian) the electrostatic energy diverges and only the 
expression (|26l) for the metric function remains valid. 
For NEDs with IRD solutions in flat space the G-ESS 
solutions are asymptotically flat, but not Schwarzschild- 
like (asymptotically anomalous behaviour) and the ADM 
mass cannot be defined [TlJ. In these cases the metric 
function takes the form (|27p if the central-field behaviour 
belongs to cases Al or A2 (but now the interpretation 
(|28p of the constant U makes no sense). For the families 
with UVD-IRD behaviour, only the expression (1251) de- 
scribes properly the form of the metrics, which are ther- 
modynamically meaningless. 

At finite r the geometric structure of the G-ESS so- 
lutions is characterized by the central behaviour of the 
electric field in flat space. The event horizons (if any) are 
located on the largest solution of the equation g(r) = 0. 
For models with UVD central-behaviour solutions the 
metrics exhibit a structure similar to that of the Reissner- 
Nordstrom solution (naked singularities, extreme black 
holes or two- horizons black holes). If the central be- 
haviour belongs to cases Al or A2 we have, in addi- 
tion, new kinds of solutions: non-extreme single-horizon 
black holes (for families Al and A2) and extreme and 
non-extreme black points (for families A2) [lOj . 

All the G-ESS solutions of the Einstein equations min- 
imally coupled to admissible NEDs must belong to one 
of the twelve classes obtained as combinations of these 
asymptotic and central-field behaviours. Consequently, 
these results, which are summarized in table HI charac- 
terize qualitatively the possible geometric structures of 
the G-ESS solutions of all admissible G-NEDs. 



IV. THERMODYNAMIC ANALYSIS 

Let us consider now the thermodynamic properties of 
the G-ESS black hole solutions of the admissible NEDs. 
In Ref.[16| D. A. Rasheed established the validity of 
the zeroth and first laws of black hole thermodynamics 
for BH solutions of large classes of field theories min- 
imally coupled to gravity, including the G-NEDs with 
Maxwellian weak field limit. In the present context the 
rotational degrees of freedom are absent and the Rasheed 
proof is restricted to the nonrotating ESS black hole so- 
lutions corresponding to our asymptotically coulombian 
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Table I. Possible geometric structures of the G-ESS solutions 
associated to the different families of admissible NEDs ob- 
tained from the combinations of asymptotic behaviours Bl, 
B2, B3 and IRD with the central- field behaviours Al, A2 
and UVD. Labels: FEFS(= finite-energy flat-space solu- 
tions), DEFS(= divergent-energy flat space solutions), AS(= 
asymptotically Schwarzschild), AA(= asymptotically anoma- 
lous), NS(= naked singularities), EBH(= extreme black 
hole), lhBH(= single-horizon black holes), 2hBH(= two- 
horizons black holes), BP(= black points), EBP(= extreme 
black points). 





Bl, B2, B3 


IRD 


Al 


FEFS 
AS 

NS, EBH, 2hBH, lhBH 


DEFS 
AA 

NS, EBH, 2hBH, lhBH 


A2 


FEFS 
AS 

NS, EBH, 2hBH,lhBH, 
BP, EBP 


DEFS 
AA 

NS, EBH, 2hBH, lhBH, 
BP, EBP 


UVD 


DEFS 
AS 

NS, EBH, 2hBH 


DEFS 
AA 

NS, EBH, 2hBH 



cases B2. As we shall see, it is easy to extend the va- 
lidity of this first law for nonrotating solutions including 
the cases Bl and B3 (the zcroth law is trivially satisfied 
for the ESS black hole solutions). However, as already 
mentioned, in the asymptotically IRD cases, for which 
the ADM mass is not defined and (as we shall see) the 
electrostatic potential diverges at infinity, the thermo- 
dynamic relations cannot be established, at least in the 
usual way. Consequently, the following considerations 
concern the B cases only. 

The set of G-ESS (nonrotating) black hole solutions as- 
sociated to a given G-NED can be assimilated to "states" 
of a system characterized by two independent parame- 
ters, the most immediate ones being the integration con- 
stants: the mass M, well defined for asymptotically-B 
NEDs, and the charge Q. The location of the black 
hole event horizons is the largest solution rh(M,Q) of 
the equation 



been extensively discussed in Ref.[lO|. The behaviour of 
M as a function of for several representative values 
of Q, obtained from Eq. (f3"0)) . is plotted in the figures [21 
and[3]for the A2 and Al-UVD cases, respectively. The 
location of the minima in these curves (when present) 
correspond to the horizon radii (fhe) and masses (M e ) 
of extreme black holes, determined by equation (|30[) and 
the extremum condition 



M 




Figure 2. M — diagram for the BH solutions of G-NEDs 
belonging to A2 families for several values of Q. There are 
three subcases, corresponding to Q = Q c = (167ra) _1 . In case 
A2a the curves exhibit minima, corresponding to the EBH 
configurations. The dotted line linking these configurations, 
together with the short-dashed straight line of Schwarzschild 
BHs (Q — 0), define the thermodynamically pertinent do- 
main, corresponding to the external BH regions. The domain 
over the former curve corresponds to the BH interiors, where 
the dashed pieces of the constant-Q lines give the M — 
relations for the inner horizons. The continuous pieces of the 
constant-Q lines correspond to the outer (event) horizons, for 
which M is always a monotonically increasing function of 
for any admissible model. When Q — » Q c (case A2c) the 
EBH configurations converge towards a extreme BP, unique 
for a given A2 NED. In the A2b case, for fixed Q < Q c , we 
can have single-horizon BHs (with M > ei) and non-extreme 
BPs (with M — ei). This picture has been obtained from 
the Bl model, as representative of the behaviour of the BHs 
belonging to the A2 families (with a > 2 in this case). 



M 



2 



(r h ,Q), 



(30) 



obtained from the condition g(rh, Q, M) = on Eq. ([26|) . 
Owing to the admissibility requirements, e ex (r^Q) is 
a monotonically decreasing and concave function of r 
at constant Q, vanishing asymptotically as r q+1 with 
q < — 1 in the B cases [8|. It exhibits a vertical asymp- 
tote at r — in case UVD or takes, in cases Al and 
A2, a finite value there (= Q 3 / 2 e(Q = 1); see Eq.flU])) 
with divergent slope in case Al. In case A2 we distin- 
guish three qualitatively different subcases, according to 
the values of the charge: case A2a if Q > Q c , case A2b 
if Q < Q c and case A2c when Q — Q c . The geomet- 
rical structures associated to these three subcases have 



dM 
dr h 



47r^T u (r ftl Q)=0. 



(31) 



The cut points between the constant-mass horizontal 
lines and the positive-slope branches of the curves of 
constant charge define the event horizon radii and the 
masses of the different black hole configurations. For a 
fixed Q, if M < M e there are no horizons and the solu- 
tions are naked singularities. Otherwise, besides extreme 
BHs, we have two-horizon (Cauchy and event) BHs, sin- 
gle horizon BHs (in cases Al and A2) and extreme (in 
case A2c) and non-extreme (in case A2b) black points. 
Moreover, since the radii of the external (event) horizons 
of the possible black hole configurations lie always on the 
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Figure 3. M — rn curves at constant Q for the Al (main 
frame) and UVD (small frame) families. The solid, dashed 
and dotted lines have the same meaning as in figure [2] The 
curves exhibit vertical asymptote at th = in the UVD cases 
and are finite there, but with divergent slope, in the Al cases. 
For the UVD and Al cases we are plotting the RN and EH 
solutions as representative examples. 



positive-slope region of the constant-Q curves, the mass- 
radius relation is a monotonically increasing function for 
all the black hole solutions of the admissible G-NEDs. It 
approaches the Schwarzschild relation M = r/,/2 at large 
Th,-, as well as at vanishing Q. 

Before going forward in this way let us define other 
important "state functions" and establish the validity of 
the first law for the ESS black hole solutions associated 
to the set of all admissible NEDs of asymptotic types B. 

From its general definition [32| the surface gravity (k) 
for these charged BHs takes the form 



l~8mjT°(r h ,Q) 
2r h 



(32) 



where Eqs. ([Tf5|) . (TT5|) and ([2TJ)) have been used. If we look 
now for the total differential of the "energy" M as a func- 
tion of Th and Q in Eq. (j30|) . and using again Eqs. lfToT) and 
(|18[) we obtain for the partial derivatives the expressions 



and 

dM 
~dQ 



dM 
dr h 



de. 



l-8wrlT°(r h ,Q) 



(33) 



dQ 



= s ; .- / E(R, Q)dR = 8wA (r h ,Q), 

(34) 

where Ao(rh, Q) is the covariant time-like component (in 
the coordinate system (f23|) ) of the four- vector potential 
in the appropriate gauge, which is fixed through the con- 
ditions 



E(r,Q) 



dr 



A (r,Q) ; A (oo,Q) = 0. 



(35) 



Thus, the total differential of M reads 



dM = K (r h , Q)r h dr h + 8irA (r h ,Q)dQ. (36) 

This is the differential form of the first law, valid for the 
black hole solutions of any G-NED of the nine families 
defined by the combinations of the three B asymptotic 
behaviours and the three central- field behaviours UVD, 
Al and A2. With the usual definitions of the "temper- 
ature" T and the "entropy" S in terms of the surface 
gravity k and the horizon area A as 



T = 7T- 



1 - 8nr*T§ 



q A 



4irr h 

together with the usual notation 



8TrA (r,Q) = $(r,Q) 



(37) 



(38) 



for the electric potential, the differential expression of the 
first law takes the more familiar form 



dM 



T{r h ,Q)dS- 
T(S,Q)dSl 



-$(r h ,Q)dQ 
$(S,Q)dQ. 



(39) 



Let us now analyze the important case of the extreme 
BHs. The condition (|31|) does not involve the mass and 
then, defines the radii of the extreme BH configurations 
as a function of the charge {rhe{Q)) which is monoton- 
ically increasing. Indeed, the general expression of the 
derivative of this function can be obtained from (|31j) and 
the first integral (|6|) in terms of the values of the field on 
the extreme BH horizon. It reads 



E, 



dr he _ 

dQ ~ rhe <p(E*y 



(40) 



and is always positive for admissible models (when Q > 
0). The explicit form of rh e {Q) which, as we shall see at 
once, determines the whole thermodynamic behaviour of 
the subclass of extreme BHs, depends on the particular 
model and can be directly obtained from Eqs.©, ([To]) 
and (|3"Tj) . once the form of the Lagrangian is specified. 
Alternatively, it can be obtained by integrating Eq.(|40|) 
with the appropriate boundary conditions, In fact this 
integration can be formally performed using the scale 
properties ([7]) leading, after some manipulations, to the 
general expression 



Q = Q ex.p[-fl(R he )} ; (R he 



The 



) (41) 



where the function VL is defined by the integral 



ft(x) 



zdz 



2E(z,l) 
•fi(EHz.l)) 



(42) 



9 



Qo and xq being related integration constants. For the 
Einstein-Maxwell theory this relation reduces to the sim- 
ple and well known expression for Reissner-Nordstrom 
extreme BHs. But for general NEDs it can become 
largely involved. Nevertheless, the The ~ Q relation for 
small and large radii is typical of each family. In both 
limits it can be determined by replacing the expression 
(TTTj) (with i = 1, 2) in ([3"Tj) and using the central or large-r 
behaviours of the field (E(r — > 0) ~ r p ; E(r oo) ~ r q ). 
This calculation leads to relations of the generic form 



The 



A; 



(43) 



where 7 j, (i = 1, 2) are the exponents defined in Eqs.([T 
or (fT3|) and the Xi are given by 



Ai = 



[87r(2 74 - 1)] 



(27i"l) 



27* 



(44) 



the ai being the coefficients in Eq. fTTj) . 

In the A2 cases the small-radius extreme BHs are 
present for charges slightly over the critical one {Q > Q c , 
see Fig|2]) and the charge-radius relation in this region 
can be obtained by a similar procedure from Eqs.® and 
ID, leading to 



The 



(2 - a)a (Q - Q c ) 



2b Q(!- CT / 2 ) 



(45) 



if a < 2, 



I 2a 

ip{a 2 ) 



(Q-Qc), 



(46) 



be immediately obtained the ones from the others. The 
following discussion about the relations between the dif- 
ferent state variables ( "equations of state" ) results from a 
(sometimes involved, but straightforward anyhow) analy- 
sis of their scale properties, their central and asymptotic 
behaviours and the admissibility conditions. 



A. The field E(r h ,Q) 

The field on the horizon as a function of r^, at constant 
Q, is obtained from the first integral ([5]). The behaviour 
of this function around r/j = and as 00 has been 

completely determined in section II for all classes of ad- 
missible NEDs and is plotted qualitatively in figure @] It 
should be stressed that the only requirement to be sat- 
isfied by the field as function of r at constant Q, in all 
admissible cases, is to be monotonically decreasing, but 
its second derivative can change sign many times at fi- 
nite r, leading to less smooth lines than the dashed ones 
plotted in this figure (see the upper small frame in figure 
IT21 below) . In determining the behaviour of E(rh,Q) as 
a function of Q at fixed we can use the scaling law (J7J . 
This leads to the condition 

E(r h ,Q) 



(r h ,Q) 




A2 



31,B2,B3 



if a > 2 and 



rl e ]n(r he ) ~ — (Q c - Q), 
oo 



(47) 



when (7 = 2. Once this The ~ Q relation is determined, 
Eq. (fBTJf gives the mass of the extreme BHs as a function 
of their radius or their charge. The same will arise for 
the other state functions defined so far, fully determining 
the thermodynamic properties of the set of extreme BH 
configurations associated to any admissible NED. 

We shall return at once to these extreme BH configu- 
rations, but let us before analyze the behaviours of the 
state variables for the general BH solutions of the dif- 
ferent classes of admissible models, as functions of r% 
(or equivalently S) and Q. This choice of the indepen- 
dent variables is made by convenience, r% being an ar- 
gument of most of the functions explicitly obtained from 
the field equations. Moreover, in the pertinent ranges 
is a monotonically increasing function of M (at fixed Q) 
and the behaviours in terms of (M, Q) and (r^yQ) can 



Figure 4. Qualitative behaviour of E(rh,Q) as a function of 
th at fixed Q. At the center, in cases UVD and Al the field 
diverges, while it takes a finite value in case A2. In this lat- 
ter case the slope of E(r, Q) at vanishing depends on the 
value of a (see Eq.((9|). The small frame exhibits the corre- 
sponding behaviours of the potential &(rh,Q)- The dashed 
pieces of the lines, matching the central and asymptotic field 
behaviours, are always monotonically decreasing but can ex- 
hibit more involved forms than those showed here for many 
models (see FigfTJ] below). The thermodynamically pertinent 
parts of these curves start at the EBH radius r^ e , obtained 
by solving Eq. (|41[) for the corresponding value of Q. 



dE 
dQ 



R h dE(R h ,Q = 1) 



2Q 



dR h 



>0, 



(48) 



(R h = -^) and we see that E, at constant 77,, is an 
increasing function of Q. Equations ([7} and (T4"8"|) allow to 
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perform a complete analysis of the behaviour of E around 
the limit Q = (ify, — > oo) and as Q — > oo (R h — » 0) 
for the different families of NEDs considered here. This 
behaviour is plotted in figure [5j 



as Q — > Q c 



Q) 

r h =0(A2) 



E(r,Q) 



1 *A2(p=0) 
1/ ■ 



31 



/I 
/ I 
/ I 

fa / I 



,EBH(A2) 



r h =0(A2) 



;ebh(ai,uvd) 




■'Al,UVD(-2<p<0) 



BJ ,<' ^Max 
2£ ' - 



Maxwell 



Figure 5. Qualitative behaviour of E(rh,Q) as a function of 
Q at constant rh, obtained from the scaling law (0 for the 
A2-B cases (left panel) and the Al-B and UVD-B cases. 
At small Q the field vanishes with three different slopes, cor- 
responding to the three possible large-r behaviours. For large 
Q, in the UVD and Al cases (right panel) the field diverges 
as Q — > 00, and exhibits different slopes depending on the 
degree of divergence at r — (p is the exponent in the first 
of Eqs.(|8}). In the A2 cases (left panel) E is limited from 
above at a fixed value (E max — a) reached at rh — for any 
Q. For < Q < Q c the dashed line E — a corresponds to the 
set of non-extreme BPs. For Q > Q c the curve of EBHs is 
plotted here for A2 cases with a > 3. These curves define the 
upper boundary (dashed-dotted) of the thermodynamically 
pertinent region. For the Al and UVD cases this boundary 
is given by the corresponding curve of EBHs (dashed-dotted 
also). The points and branches over these boundaries corre- 
spond to the inner horizons. The matching between the small 
and large Q parts of a constant-r^ curve depends on the par- 
ticular model, but is always monotonically increasing. The 
dashed piece of the straight line in the right panel represent 
this matching for the Maxwell (RN) case. The dashed piece of 
the curve in the left panel corresponds to the similar matching 
for a particular A2-B3 model. 



The potential <&(r h ,Q) 



The expression of the electrostatic potential as a func- 
tion of rh and Q can be easily obtained from that of 
E(rh, Q) through the relations (|55|) . It is plotted at con- 
stant Q in the small frame of figure SJ In looking for the 
behaviour of $ as a function of Q at constant r^ we first 
obtain its scaling law from Eqs.([7]) and (|35|) . which reads 



(51) 



$(r h ,Q)= y/Q$(R h ,Q = l), 



(R h = -^|=) and shows that $ is a homogeneous function 

of r^ and y/Q of degree one. From this law we obtain the 
Q-derivative of $ as 



5$ 
dQ 



$(r h ,Q) + 8irr h E(r h ,Q) 
2Q 



>0, 



(52) 



which shows that $ is a monotonically increasing func- 
tion of Q at constant r/, . The form of this derivative and 
of the scaling law allows straightforwardly for a complete 
discussion of the behaviour of the potential at small and 
large Q for constant r^. This behaviour is displayed in 
figure [6] for the different families. In the A2 case, for 
Q < Q c , the curve r/, = corresponds to the non-extreme 
BPs and defines a part of the upper boundary of the ther- 
modynamically pertinent region. The Q = Q c point of 
this curve defines the extreme BP. For large values of Q 
the potential at constant behaves as 

$(r ft , Q -> 00) ~ y/Q${r h = 0, Q = 1) - 8vrar ft (53) 
in the A2 cases. In the Al cases we have 



For the extreme BHs we can obtain easily the be- 
haviours of the field on the horizon at small and large Q 
by replacing the r^e ~ Q dependencies obtained for the 
different families. The large Q behaviours and the small 
Q behaviours in the cases UVD and Al are obtained 
from Eg. (1431) and have the same form, which reads 



E e {Q) 



iiQ~ 



(87r)( 2 ^-i)(2 7l -l)' 



(49) 



as Q (i = 1) and as Q —> 00 (i — 2). In the same way, 
for the A2 cases the behaviour of E e (Q) for small-radius 
extreme BHs (Q — > Q c ) is straightforwardly obtained 
from equations (j44 |) -(|47 |) in models with different values 
of it. For example, for a > 2 we have 



E e (Q)~ a -b 



2a 


a/2 


'Q-Qc 


cr/2 


_(p(a 2 )_ 




[ Q 





(50) 



$(r & ,Q-»-oo)~v / Q*(»Vi = 0,Q = l) (54) 
whereas in the UVD cases we have 



$(r h ,Q -> 00) - -8tt 



p+i 



(55) 



(p < — 1). At small Q the behaviour of the constant 
rh curves is dominated by the large-r behaviour of the 
electrostatic field, and reads 



*(r h ,Q->-0) 



q+1 



(56) 



where -1 > q > -2 (Bl), q = -2 (B2) and q < -2 
(B3). 

For Q > Q c the thermodynamically pertinent region 
lies under the curve of extreme BHs, which is defined by 
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eliminating r^ e between equation (|31[) and the expres- 
sion of $(r^,Q). The shape of this curve, which begins 
on the extreme BP (Q = Q c ) in the cases A2, and at 
Q = in cases Al and UVD, depends of the particular 
model, but its behaviour around these points and also 
for large values of Q is universal for every family. It can 
be determined explicitly from equations (|43l) to (|47p and 
the corresponding expansions of $(r, Q) around r — > 
or r — > oo, obtained from the integration of Eqs.© and 
©. For large values of Q these equations lead to 



1' 



S(0,Q) = 


VcTs (o,i) 




r h = 0^ 






/ -' 


f ,--"5 

A' 
\ 


p<-2y 




\ 


P ,?'i'2<p<-l 




EBH^ 












Qc 



< —1). We see that the potentials on the horizon for 
large-Q extreme BH configurations diverge slower than 
y/Q in the Bl cases, approach a constant value y/8Tv / ' a 2 
(the same in all the asymptotically coulombian (B2) 
cases if the normalizing ai constants of the Lagrangians 
in Eq. lfTTT) are fixed to the same value) and vanish in the 
B3 cases. For small Q values in the Al cases the electro- 
static potential on the horizon of extreme BHs behaves 
as 

$ e (Q^0)~$(0,l)VQ- (59) 

In the UVD cases the same expressions ((57)1 and (|5"5)) 
hold, with the replacements q — ► p, «2 — > &i, 72 — > 7i 
and fj.2 — > Mi (excepted in cases with p = —1, where 
logarithmic dependence arises, but does no modify the 
behaviours described here). Now, $ e (Q) vanishes at Q = 
with divergent slope for —2 < p < — 1 cases, takes a 
finite value there (= y^njai) if p = —2 and exhibits a 
vertical asymptote if p < —2 (see figure [7J. 

In the A2 cases the integration of Eq. (|9|) , together with 
(|I5])-([1| lead, for a < 2, to 



Figure 6. Qualitative behaviour of &(rh, Q) as a function of Q 
at constant Th, for small and large values of Q, deduced from 
the scaling law (|51[) for the A2 cases (main frame) and for 
the UVD and Al cases (small frame). As Q — > the form of 
the constant curves is governed by Eq. H56[) . deduced from 
the asymptotic behaviour of E(rh — > oo, Q) (B-cases). The 
upper curve in the main frame ($ = V / Q$(0, 1)) gives the be- 
haviour of the potential on the horizon of the BPs (rn = 0) for 
Q < Q c in the case A2. The large-Q shapes of the constant- 
Th curves are determined by the central behaviours of the 
field E(r h — > 0, Q). In cases A2 these curves behave asymp- 
totically as parabolic branches, paralel to the upper parabola 
(see Eq.([53}). The asymptotic shapes of the Al and UVD 
curves depend on the degree of divergence of the field at the 
center (see equations (|54[1 and (|55j). The form of the dashed 
pieces of constant rn curves, matching both regimes, depends 
on the particular model but must be always monotonically in- 
creasing. In fact, the thermodynamically meaningful region, 
where the external horizons are present, lies under the = 
curve for Q < Q c and the line of EBH solutions for Q > Q c 
in the A2 cases. For the Al and UVD cases this region lies 
under the line of extreme BHs, which is calculated in the text 
for the different cases. It is qualitatively represented in the 
main frame for Q > Q c in a particular A2 model with a > 1 
(dashed-dotted line). The small-Q, Q > Q c and large-Q be- 
haviours of these EBH curves defining those upper boundaries 
are plotted in Fig[7J for all possible cases. 



$ e (Q^oo)-^ 2 Q 2 " 
where the positive coefficient is 



M2 



2-q (2- q\« 



4(l + g)Va7yJ V 327r 



(57) 



(58) 



where 



$(o,i)- w (i-^ Nl/ ° 



, (60) 



oj = ana 



l + l/a f 2 -< J 

2b 



l/a 



(61) 



whereas for a > 2 we obtain 



*e(Q -> Qc 



$(0,1) - 87T, 



' 2a 3 
<p(a 2 ) 



(62) 

The derivative of § e (Q) can be directly obtained from 
Eg. ([5^)1 and PU|) in the general case and reads 



dQ 



1 

2Q 



<L> r 



E„ 



r he <p(E 2 e ) 



(63) 



This equation allows the calculation of the slope of the ex- 
treme BH curves around Q c , which is positive and finite 
in the cases a < 1, finite positive, null or finite negative if 
a = 1 and negative divergent if a > 1. All these features 
are displayed in figure [7J 

The function $ e (Q) involves the two thermodynam- 
ically conjugate state variables $ and Q, and defines a 
vanishing temperature "isotherm" of the set of BH states 
associated to a given model (see equation (l37l) and the 
extremality condition (I5T|) '). Thus the analysis of these 
curves can reveal the presence of van der Waals-like first- 
order phase transitions. This requires, as a necessary 
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/Qc S(0, 1) 



2<q<-l (Bl) 




for the event horizons of the BH configurations in the 
A2 cases is bounded by above by the line of extreme BH 
configurations and the line — 0, which intersect in a 
point with Q = Q c , corresponding to the extreme black 
point configuration A2c. In case Al the line of extreme 
BHs is always below the line = and defines the upper 
boundary of the BH region. This case is obtained from 
FigjS] in the limit of divergent maximum-field strength 
(a — >• oo or Q c = 0). For the UVD case (in the small 
frame) this boundary is also given by the line of extreme 
BHs. The Q axis, corresponding to — » oo, is the lower 
boundary of this thcrmodynamically pertinent region. 



Figure 7. Behaviours of the potential on the horizon of the 
extreme BHs (<E> e (Q)) for small values of Q (in the Al and 
UVD cases), large values of Q (in the B cases) and Q > Q c 
(in the A2 cases). See the discussion in the text for details. 



condition, the existence of some values of the charge for 
which the slope of $ e (Q) vanishes. Thus, the search for 
zeroes of ([63)) is a first step in the identification of such 
transitions. A glance to figure [7] shows that these spe- 
cial values of the charge must be necessarily present for 
models belonging to families with Al or UVD central 
behaviour (with —2 < p < 0) and B3 asymptotic be- 
haviour. The same arises for models belonging to the 
UVD (with p < — 2) and Bl families, or for those be- 
longing to A2-B1 families with a > 1. Other interesting 
cases are the models belonging to the A2-B2 family with 
a > 1, when the condition 

$(0,1) < 8ttW— (64) 

V Oi2 

holds. But, beyond this necessary condition, the confir- 
mation of the presence of such phase transitions in a given 
model requires the simultaneous analysis of the Gibbs 
free energy, defined as 

G = M — - TS = Th ~ M (65) 

for each set of extreme BHs. A systematic exploration of 
these issue is in progress. 

C. The external energy function e e x(rh,Q) 

As a function of r^, for fixed Q, this function has been 
extensively discussed and employed in Ref.[l(| for the 
characterization of the different families. Its typical qual- 
itative behaviour is plotted in figure [5] for all cases. The 
dependence of e ex on Q at constant r% can be determined 
from the scaling law (|20]l and Eq. (f34"]) . This leads to the 
curves plotted in figure^ corresponding to the cases Al- 
A2 and UVD (small frame). The representative region 



£ e x (r h , Q) 
£ (Q) 



1 




Figure 8. Qualitative behaviour of e ex (r, Q) as a function of 
Th at constant Q for the Al, A2 and UVD cases (in units 
of the total energy e(Q) in the Al and A2 cases). The three 
A2 curves correspond to three values of the charge Q > Q c 
(cases A2a), Q < Q c (cases A2b) and Q = Q c (case A2c). 
The cut points of the (dashed) straight lines defined by y = 
1 — rh/(2e(Q)) with the A2 curves give the location of the 
horizons (see Eq. (J30J) ) ■ 



D. The Mass M(r h , Q) 

As a function of at fixed Q the behaviour of the 
mass has been already discussed at the beginning of this 
section and plotted in Figs. [5] and [3] The behaviour of 
M as a function of the charge at fixed rh can be now 
obtained from the scaling law 

M(r h , Q) = (1 ~ f^ Rh + Q"/ 2 M(R h ,Q = 1) (66) 

[R h = which results from (|30|) and the last of 

Eqs. ([20l) . This function is plotted in figure [TO] for the 
cases A1-A2 and UVD (small frame). The curve rep- 
resenting the extreme BH states is also displayed on this 
figure. It is the envelope of the beam of constant rh 
curves (see Eq. (|3T|) ) and begins tangent to the = 
curve at Q = Q c in case A2. In case Al the correspond- 
ing curve is obtained by moving this tangency point to 
the origin (Q c = 0, in the limit a — > oo). In cases Al 
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Figure 9. Behaviour of e e x(rh,Q) as a function of Q at con- 
stant rh for the A2 (main frame) and UVD (small frame) 
cases. The sets of EBH configurations are drawn in both 
cases. The beams of parabolic-like curves correspond to dif- 
ferent constant values of th- The points of the region under 
the composite curve = and EBH in the case A2, and 
under the EBH curve in the UVD case correspond to the 
exterior of the BH configurations, whereas the points of the 
regions over these curves are associated to the BH interiors, 
playing no role in the thermodynamical analysis. The typical 
behaviours in the Al cases can be obtained from those of the 
main frame figure by allowing the cut point of the ru = and 
EBH curves to go to the origin, i.e. by making Q c — > 0. In ob- 
taining these curves we have used the BI and the RN models 
as representatives of the A2 and UVD cases, respectively. 



and A2, the region over the = curve contains the 
single-horizon BH states, the region between this curve 
and the extreme BH curve contains the two-horizon BH 
states and the region under these curves and over the 
Q axis corresponds to naked singularity configurations. 
Similar comments can be done for the UVD case (small 
frame in figure I10p , where the beam of constant-radius 
curves approach the M axis as Th — > 0. 



E. The temperature T(rh,Q) 

In analyzing the behaviour of the temperature Q37p as a 
function of at constant Q let us introduce the function 



v(rh,Q) 



1 dM 



1 -87rrgr° (r h ,Q) 

47T 



= r h T(r hl Q), 
(67) 



2tt dr h 
or, equivalently, 

l(Th,Q) = ^ - ±QE(r h , Q) + 2r 2 h cp(E 2 (r h , Q)), (68) 



where the last hand side of Eq. ([T6f has been used. The 
derivative of this function is easily obtained and, as a 
consequence of the positivity of (p, satisfies 



Figure 10. Behaviour of M as a function of Q at fixed 
fh for the A2 (main frame, BI model) and UVD (small 
frame, RN solution) cases. Some constant ru curves are 
shown and are split (excepting for the rn = one) into a 
continuous part, corresponding to the event horizons, and a 
dashed part, corresponding to the inner horizons. The EBH 
configurations are the envelopes of the beam of constant 
curves. The = curve in the A2 cases (corresponding to 
M = e(Q) = Q 3/2 s(Q = 1)) defines the non-extreme BPs 
for charges below Q c and the EBP configuration for Q — Q c . 
For Q > Q c this curve separates the regions associated to 
the single horizon and two-horizons BH configurations. The 
diagram corresponding to the Al cases is similar but with 
vanishing critical charge (Q c = 0). The analysis of the dia- 
gram in the small frame (UVD cases) is similar, but now the 
th — > curves approach the Q = axis (which corresponds 
to the Schwarzschild BH limit) and the charged single-horizon 
BH configurations are absent. 



dr] 
drh 



1 d 2 M 



2n drl 



4r h ip > 



(69) 



Thus 77, at fixed Q, is a monotonically increasing func- 
tion of rh for any admissible G-NED and, owing to the 
behaviour of ip and E as, Th — > 00, approaches asymp- 
totically the constant 1/47T. The behaviour of this func- 
tion is plotted in figure [TT] for representative examples 
of cases Al (Euler-Heisenberg) and UVD (Reissner- 
Nordstrom) for a given value of Q. The parts of these 
curves in the negative rj region correspond to the inner 
horizons and must be discarded for the thermodynami- 
cal analysis. The straight lines drawn from the origin to 
the points of these curves in the positive 77 region have 
slopes tan(#) = T, which define the corresponding tem- 
peratures. Since these curves are continuous, monoton- 
ically increasing and approach asymptotically the con- 
stant value rj — 1/47T for any Q, there is always an abso- 
lute maximum value of the angle 9. Thus, the tempera- 
ture ranges from zero, on the cut points with the Th axis 
(corresponding to the extreme BH configurations), to a 
maximum value T max , corresponding to an angle 9 max 
in Fig llll This maximum (and, eventually, other local 
extrema of the curve) correspond to BH configurations 
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whose horizon radii are solutions of the equation 



4 71 




us illustrate this behaviour with a particular example. 
Consider a family of NEDs whose ESS solutions are of 
the form 



(73) 



where R — r/^/Q, the parameter Q being the charge 
associated to a given solution and a > 0, /3 > and the 
integer n > 2 being fixed constants for a given model 
of the family. These fields are monotonically decreasing 
functions of r if the parameters fulfill the (sufficient but 
not necessary) condition 



Figure 11. Behaviour of rj(rh,Q) as a function of th and 
a fixed value of Q for the UVD (RN) and Al (EH) cases. 
The slopes of the dashed straight lines matching the origin 
with the points of the r\ curves give the temperatures of the 
associated BH configurations. The r\ = points correspond 
to vanishing temperature EBH configurations. The associated 
temperature behaviours are plotted in the small frame. 



dT 
dr h 



4r%ip - i] 



T 

Au> = 0. 



(70) 



which results from Eqs. (|57|) - (IM|) . Moreover, from 
Eq. (1671) and the asymptotic behaviour of the energy den- 
sity in the B-cases it is easily seen that the temperature 
vanishes as 



T(r h -t oo, Q) 



0. 



(71) 



for any large mass BH (this behaviour is related to the 
asymptotically Schwarzschild character of these space- 
times). The BH temperature as a function of at fixed 
Q is plotted in the small frame of figure QTJ The slope of 
this curve is given by the derivative defined in Eg. ([70]). 
The temperature as a function of M has a similar shape, 
but the slope of the T — M curve, obtained from Eqs. (f6"Tj) 
and (ITUI). reads 



dT 
dM 



Q 



4r h ip - T 
2<Kr\T ' 



(72) 



and diverges for the extreme BH configurations. Owing 
to the convex character of the curves t](rh, Q) in the par- 
ticular cases of figure [TT] there are at most two cut points 
with the radial lines for T < T max or, equivalently, there 
are at most two roots of Eq.(|70p in these cases. How- 
ever, in more general cases, the sign of the curvature of 
vfrh, Q) a t constant Q can change many times as in- 
creases, leading to many cut points with the radial lines 
or, equivalently, to many solutions of equation (|70|) . Let 



a < 



6 



(74) 



and correspond to the ESS solutions of charge Q associ- 
ated to a family of admissible NEDs whose Lagrangian 
densities are parameterized by the constants a, /3 and 
n. The form of these Lagrangian densities can be di- 
rectly obtained from Eqs.Q and (|73j) through a quadra- 
ture. Owing to the central and asymptotic behaviour 
of the fields, the family belongs to the UVD-B2 cases 
(as the Maxwell Lagrangian). The form of this field is 
plotted in the upper small frame of FigfTS] for the model 
with the particular choice of the parameters n = 4 and 
a = \/f3= 1/5.2, satisfying the condition (|74p. and for 
the characteristic value of the charge Q = 1. In the 
same figure we have plotted the corresponding function 
r l{ r h, Q) and the associated temperature T{rh, Q) (lower 
small frame). As can be seen there are in this example 
three solutions of Eq. (f70j) corresponding to a local maxi- 
mum, a local minimum and an absolute maximum in the 
T — rh diagram. Because of the monotonic relation be- 
tween rh and M, a similar behaviour must be obtained 
in a T — M diagram, where the slope of the temperature 
decreases from +oo at the extreme BH configuration and 
vanishes on the first local maximum, defining a stable 
phase in this range. In the intermediate region between 
the local minimum and the absolute maximum, this slope 
becomes again positive, defining a second thermodynam- 
ically stable phase there. For increasing values of n in the 
family (fT3"|) the number of solutions of Eq. flTO")) increases, 
and new ranges of between extrema support stable 
phases. We conclude that for models whose ESS black 
hole solutions belong to cases UVD and Al (with any 
B-type asymptotic behaviour) the temperature increases 
always from zero, for the extreme BH configurations, and 
vanish asymptotically for large mass BH configurations, 
as for the RN black holes. In the range of intermediate 
masses the temperature has always an absolute maxi- 
mum and can exhibit a set of local extrema, determining 
a more or less involved phase structure, which is fully 
characterized by the solutions of Eq. ([70|) . 

The behaviour of the temperature as a function of the 
charge at constant can be analyzed from its scaling 
law 
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4717] (r h ,Q) 




Figure 12. Behaviour of the curve r/(rh,Q) as a function of 
fh for BHs obtained from the gravitating field of Eq. (|73|) . 
whose profile is shown in the small upper frame. Note that 
the change in the sign in the curvature of the field leads to 
the existence of up to three radial straight lines tangent to 
the curve ri(rh,Q), which define an absolute maximum and 
two local extrema of the associated temperature curve (small 
lower frame). 



. 4 



0.2 




Figure 13. Behaviour of T(rh,Q) as a function of Q at con- 
stant m for the Born-Infeld model, as representative of the 
A2 cases. The set of EBH states lies on the Q axis for Q > Q c 



To address this question let us go to the next order in 
the polynomial expansion (J9j) of the field solutions of a 
given A2 model around the center, which can be written 
as 



T(r h ,Q) = 



1 - 



y/QT(R h ,Q = l) 



(75) 



(R h = -^) which can be directly deduced from Eq. ([rJ5|) 

and the first integral ©. From this equation we obtain 
the derivative 



dT 
dQ 



(76) 



which is negative for any Q > and vanishes at Q = 0, 
where the temperature reaches the Schwarzschild value 
T = (Anrh)" 1 - Figure ITUl displays the temperature as a 
function of the charge at different values for the case 
of field solutions E(r) with monotonic slope (Born-Infeld 
here, as representative of A2-B2 cases). 



F. The cases A2 and the black points 

Although the validity of the analytic discussion of the 
precedent subsection is general, it has been applied there 
mostly to the thermal behaviour of BH solutions of UVD 
and Al families only. Let us discuss here the cases A2, 
which lead to more involved (and richer) thermodynamic 
structures. The method used in determining the be- 
haviour of the temperature as a function of the horizon 
radius (or of the mass) at constant Q still holds, but in 
this case the functions rj(rh,Q) exhibit new characteris- 
tics which lead to configurations with qualitatively new 
thermodynamic properties. The main difference with the 
previous cases concerns the behaviour of these functions 
as Th — ¥ 0, which defines the black point configurations. 



E{r) ~ a - b{Q)r a + c{Q)r x =a- b R a + c Q R x , (77) 

(R = -^). Here A > a and Q x/2 c(Q) = c(Q = 1) = Co is 

a universal constant of the model. From Eq. (|68p we see 
that the behaviour of r)(rh ~ 0, Q) can be obtained from 
Eq.flZZJ) and the behaviour of ip{r h ,Q) = ip(E 2 (r h ,Q)) 
near the black point configurations (r/, ~ 0). Using ([77)1 
and the first integral © we are led to 



V(r h ~ 0, Q) ~ 2Q 



Ac(Q)x_ a «r6(Q) __ s 



A- 2 h a-2 h 
if a, A 7^ 2. When A = 2, <p behaves as 



<p(r h ~0,Q)~ 2Q 
and if a = 2 

<p(r h ~0,Q) ~ 2Q 



2c{Q)]n{r h ) 



ab{Q) 
a-2'' 



.(7-2 



A, 
(78) 

-A, 
(79) 



HQ) x_s 
A-2 7 '' 1 



26(Q) \n(r h ) 



+ A. 

(80) 

In these equations the constants A are independent of 
the particular solution. If a > 2 they are given by 
A = ip (a 2 ) > 0. Otherwise they can be explicitly cal- 
culated from the form of the Lagrangians, as functions of 
the exponents a < 2, A and the associated universal co- 
efficients in the polynomial expansion of the field around 
the center (this may require to go beyond the order of 
the expansion in Eg. (1771) ). The corresponding expression 
for the A2 cases with a < 2 and A > 2 is 
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A = lim 

X-+a 2 



X)" 



(81) 



while for a = 2 we have 



lim 



<p(X) + 2b ln(a - VX) 



(82) 



The expression of A becomes more involved when A and 
some of the next-order exponents in the polynomial ex- 
pansion of E(r ~ 0) are smaller than two, but its calcu- 
lation remains straightforward. 

Equations (j77|) - (l82j) replaced in (|68| give the form of r\ 
for the small radius black holes in the different A2 cases. 
For the black points (jh = 0), rj takes always the finite 
value 



r){r h =0,Q)=4a(Q c -Q) 



(83) 



which vanishes for extreme BPs (Q — Q c ) and is pos- 
itive for non-extreme BPs (Q < Q c ). If Q > Q c , the 
small values of r?i(< The) correspond to inner horizons, 
for which 77(0, Q) < does not define a meaningful tem- 
perature (see Fig 03}. Obviously, for non-extreme BPs 
the temperature diverges in all cases, but the tempera- 
ture of the extreme BPs is now dependent on the value of 
the parameter a of each A2 model. Indeed, from the ex- 
pansions of rj at small rh in the different cases we are led 
to the several possible expressions for the temperature of 
the small-radius BHs: 
If <t, A 7^ 2 we have 



T(r h -> 0, Q) ~ fe±M) „ 1 - Q/Qc _ W) rf -i 



+ 



Th 47rr/ l (7-2 h 

HQ) A _ x 



A-2 



-r^ L + 2Ar h . 



(84) 



If a = 2, we have 



rpi _^ an , v(r h ^0,Q) l-Q/Qc , 8Qc(Q) x _, 
1 ' l '^^ r, ~ 4nr h + A - 2 ^ 

+ ±Qb(Q)r h {l - 2 ln(r fc )) + 2Ar h , (85) 



and if A = 2, 



T(r h ->0,Q) 



y(r h ^0,Q) ^ 1 - g/g e _ 8Qb(Q) ^_ x 

r h ~ 4 7 rr /l £7 - 2 ^ 

4Qc(Q)r h (l-21n(r h )) + 2Ar h . (86) 



In all these expressions the first term dominates if Q ^ 
Q c , in agreement with Eq. (|83|) . For the extreme BPs the 
temperature is given by the slope of rj(rh, Q c ) m r h = 
and is null if > 1, finite if £7 = 1 or divergent if £7 < 1 
(see FigE). 



4JTr7(r h ,Q) 
1 



A2b(Q<Qc) 
A2c(Q = Q c )- 
A2a(Q>Q c ) 



f 0=1 

^ a>l 




Figure 14. Small and large rh behaviour of r/(rh,Q) at con- 
stant Q for the A2 cases. The three beams formed by con- 
tinuous pieces of curves at small rh correspond to the three 
subcases defined by the values of the charge Q ^ Q c . For 
each beam the three curves correspond to different values of 
the parameter a ^ 1, exhibiting different slopes (vanishing, 
finite, or divergent, respectively). The asymptotic behaviour 
of these curves is obtained from Eqs. (|68|) . (© and (|11|) and 
is also displayed with pieces of continuous lines at large rh- 
As in the previous cases (Al and UVD), the form of the 
dashed lines (intermediate region) depends on the particular 
model and can exhibit more involved shapes than the ones 
displayed in this figure. The small figure exhibits the form of 
the temperature for the BI model (for which a = 4) for several 
values of Q. Note that the temperature of the EBPs (cases 
A2c, with Q = Q c ) is given by the slope of the corresponding 
curve at rh = and vanishes in models with a > 1, takes a 
finite value in models with a = 1 and diverges in models with 
a < 1. The temperature of the non-extreme BPs (A2b cases, 
for which Q < Q c ) is always divergent. 



It is now evident that the temperature behaviour in 
case A2a (77(0, Q > Q c ) < 0) is similar to the ones of 
the cases UVD and Al (see the small frame of Fig fT^)) . 
exhibiting an extreme BH configuration with vanishing 
temperature at finite rh, vanishing temperature BHs for 
large rh, an absolute maximum temperature BH and (de- 
pending on the particular model) many possible BH con- 
figurations corresponding to local extrema of T for some 
values of r/, (or M), which define stable phases in some 
intermediate ranges. This behaviour is qualitatively in- 
dependent of the values of the parameter £r. 

In case A2b (77(0, Q < Q c ) > 0) there is a continu- 
ous spectrum of single- horizon BHs, whose radii increase 
from r h = (when M B p = e(Q) = Q 3/2 e{Q = 1), cor- 
responding to a non-extreme black point; see figure [2]) 
to rh — > 00 as M — > 00. In these cases the temperature 
diverges for the black points and vanishes asymptotically 
for large mass BHs. In the intermediate region, depend- 
ing on the structure of the ESS fields, the temperature 
can exhibit local extrema defining ranges of stability, as 
in the example (|73p . Again, this behaviour is qualita- 
tively independent on the values of the parameter a. 

In case A2c we have also a BH spectrum going from 
rh = (corresponding now to an extreme black point) 
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to Vh — > oo. The mass of the extreme BPs equals the 
soliton energy and can be written as 



M BPextr = Ql /2 e{Q = 1 



$(r = 0,Q = 1) 
96(vra) 3 / 2 : 



(87) 



where the last term is obtained from Eqs. (f21j) . ()35|) and 
(|38|) . and is a universal relation for A2 models. The 
thermal behaviour of the extreme BPs for the different 
models is characterized by the values of the parameters 
a and A and can be directly obtained from Eqs. (lMl) - (fc>o]) 
and their derivatives: 

• For models with a < 1 the temperature diverges for 
the extreme BPs, exhibiting a vertical asymptote at 
r h = 0. 

• For models with a = 1 the temperature of the ex- 
treme BP takes the finite value (see Eq.(| 



Table II. Thermodynamic behaviour for extreme BPs (th = 
and Q = Q c ). See Eqs. (|77) l-l |82) l for the meanings of the 
parameters. 





a < 1 


a = 1 


1 < a < 2 


<T > 2 


T(r h = 0) 


+00 


V7TO 








dT 1 


— oo 


-oo(c > 0, A < 2) 
+oo(c < 0, A < 2) 
2A(A > 2) 


+oo 


2A 


S 2 T 1 

h 


+oo 


co ■ oo(A < 3) 
64co y / na(X = 3) 
0(A > 3) 


— oo 


-oo(2 < a < 3) 
— 646o\/'™'(' :r = 3) 
0"(cr > 3) 


Cq, r h =0 


0" 


0"(A < 3, A = 0) 
1K 6 ° (A-3,A-0) 
+oo(A > 3, A = 0) 
0- A 


0+ 


0+ 



T(r h = 0,Q c ) 



cHQc) 



2b 



which is a universal constant for each A2 model. 
The slope of the temperature curve at vanishing Th 
diverges for A < 2, being =Foo for cq ^ 0, respec- 
tively. For A > 2 this slope is 2A, which is finite. 
Here the sign of A determines the stability of the 
BPs. 

• For models with 1 < a < 2 the temperature of 
the BPs vanishes and its slope at = is always 
positively divergent. 

• For models with a > 2 the black point temper- 
atures vanish and the slope of the temperature 
curves at the origin is 2A > 0. 

A similar analysis of the second derivatives of the tem- 
perature with respect to Th, obtained from equations 
(|M| - (|86[) . give us the concave or convex character of the 
T — Th curves for configurations around the BPs. A sum- 
mary of these results for the extreme BP solutions in this 
A2c cases is given in table [TTJ 



G. The specific heats Cq(rh,Q) 

The specific heat at constant charge of the BH config- 
urations is defined as 



there vertical asymptote. From Eqs. (|67 
obtain the general formula 



C Q (r h ,Q) 



Ar h ip - T" 



and (l68t we 



(90) 



and eliminating the Lagrangian function ip in terms of 
the temperature and the field we are led to 



C Q (r h ,Q) 



4ir 2 rlT 



2-KTh.T + WnQE - 1 



(91) 



We see from this last expression and the asymptotic be- 
haviour of T, given by Eq. (fTTj) . that the specific heat of 
large mass BH solutions of any admissible B-type model 
behaves as 



CQ{r h -)■ oo,Q) 



-oo, 



(92) 



as should be expected from their asymptotically 
Schwarzschild character. The slope of Cq as a function 
of rh reads 



dC Q 



dr h 



Q 



2nrh 



4irr h T 

Ar h (p - T 



4Q|f+r/2 
4r h (p - T 



(93) 



C Q (r h ,Q) 



dM 



dM 
9r h 


Q _ 


2-KJ] 


dT 

9r h 


Q 


dT 1 
dr h \Q 



(89) 



and we see that, as a function of or M, diverges on the 
relative extrema of the temperature function, exhibiting 



Using the asymptotic behaviours of the different func- 
tions in this formula for the B-cases we see that for large- 
r/i BHs this derivative diverges as — 47rr^, consistently 
with the asymptotic behaviour of Cq in Eq. (|92|) . 

For the extreme BHs associated to the cases UVD, Al 
and A2a we have rj = T = and | = 4tp > (see 

Eq. flSi?)) ). Consequently CQ- ex t r = for any extreme BH 
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configuration of admissible models. The slope of the spe- 
cific heat as a function of , given by Eq. , converges 
to the value 



For extreme BPs in models with 1 < a < 2 the 
specific heats vanish, behaving as 



dC Q 



dr h 



(94) 



for the extreme BHs of any admissible model. The slope 
of the specific heat as a function of the mass is obtained 
from Eq.(|93]l as 



dC Q 



dM 



dC Q 




dr h 


Q 


dM 




dr h 


Q 



1 dCr 



2nr h T dr h 



(95) 



and diverges in the extreme limit as 1/T. 

Let us analyze now the specific heats of the BP con- 
figurations discussed above. Using Eq. (|9Tj) and replacing 
there the expressions ([54]) and (1771) for the expansion of 
the temperature and the field around rj, ~ we are led 
to 



C Q (r h ^0,Q c ) 



2n 



a - 1 

for small horizon radius BHs. 



(101) 



• For extreme BPs in models with a < 1 the specific 
heats vanish, behaving as 



C Q (r h ^0,Q c ) 



2tt 



a - 1 

for small horizon radius BHs. 



(102) 



• For extreme BPs in models with a = 1 the domi- 
nant terms in the expression of the specific heat of 
small radius BH configurations reduce to 



where 



(96) 



C Q (r h ^0,Q c ) 



2irb(Q c )(X-2)r h 



c(Q c )(A-l)r£- 2 + 47ra(A-2)A' 



when A 2 and 



(103) 



V J i 1 - -If- | r 



Qc 

64tt 2 Qc{Q) 



.A+2 



Mv 2 Qb{Q) a . 



(97) 



and 



9. 

32ttQ(A 



+ 



327r(cr- l)Qb(Q) 



l)c(Q) 



A 



+ 8?rA7 



hi 



(98) 



which is valid for a, A =^ 2. It is straightforward to obtain 
two similar expressions for the cases a = 2 and A = 2 
from Eqs. (|85p and (|86]). The analysis of these equations 
leads immediately to the following conclusions: 

• The specific heats vanish for any non-extreme BP 
(A2b cases, Q < Q c ), behaving as 



C Q (r h ^0,Q) 
for small horizon radius BHs 



(99) 



• For extreme BPs (A2c cases, Q = Q c ) in models 
with a > 2 the specific heats vanish, behaving as 



C (rfc ->0,Q c )~27rr£->0+ 
for small horizon radius BHs. 



(100) 



CQ(r h -> 0,Q C 



2ixb{Q c )r h 



c(Q c )\n(r h ) V a2c \n(r h ) 



7T b r h 



-> o, 

(104) 

when A = 2. As easily seen, the specific heats of 
extreme BPs in these a = 1 cases vanish, excepted 
for models with A > 3 and A = 0, as we shall 
see at once. If A — 2, Cq behaves for small horizon 
radius BHs as in Eq. (|104|) . where the sign near zero 
is given by that of the coefficient cq. This behaviour 
is not dependent on the value of A. For A > 2 and 
A 7^ 0, the specific heat of small-r^ BHs behaves 
as 



C Q (r h ^0,Q c ) 



2b r h 
a A 



0. 



(105) 



where A determines now the sign of the approach 
to zero. For 1 < A < 2 the specific heat behaves as 



C Q (r h ^0,Q c ) 



27t6q(A-2) ™. 
co(A - 1) 



(106) 

which vanishes with the opposite sign of Co and is 
not dependent on the value of A. Finally, for mod- 
els with A = the expression (|106j) of Cq remains 
valid for any value of A 7^ 2. In these cases we can 
have extreme BP solutions with finite specific heat 
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C Q (r h = 0,Q c ) = 



b 



IGclcq ' 



(107) 



if A = 3, besides models with divergent specific heat 
extreme BPs if A > 3. 

In the last line of table |TT] we have included these re- 
sults on the possible behaviours of the specific heats for 
extreme black point solutions associated to the family of 
A2-type gravitating NEDs. This discussion exhausts the 
possible thermodynamic structures of the G-ESS black 
hole solutions of admissible NEDs minimally coupled to 
gravity. 



H. Some examples 

In order to illustrate different points of this analysis 
let us consider some particular examples belonging to 
the different families of admissible models. 

The thermodynamic properties of the Reissner- 
Nordstrom solution of the Einstein-Maxwell equations 
are representative of the typical behaviour of the small- 
radius BH solutions of the UVD family. Although this 
RN solution has been exhaustively analyzed in the liter- 
ature (see [l4[ and references therein) let us mention a 
few of its features, owing to this representative character. 
We have already make use of this example in the small 
frames of figures [3] and [IT] representing the generic be- 
haviour of the mass and temperature as functions of the 
BH radius for this kind of solutions. We have also dis- 
played the temperature and the specific heat as functions 
of the mass and of the horizon radius in figure [15] The 
explicit expressions for these variables in the RN case 
are straightforwardly obtained from the formulae of the 
precedent analysis and lead to the expected thermody- 
namical behaviour for the small and large horizon radius 
BHs of the UVD-B2 cases. 

As a representative example of the Al family let us 
consider the Euler-Heisenbcrg effective model of quantum 
electrodynamics [f| [9[ .whose gravitating counterpart has 
been analyzed in Ref . [33[ • The interest on this model lies 
in the description of quantum corrections induced on the 
electrostatic field by the vacuum and their effects on the 
gravitational field generated by point-like charges, which 
is classically described by the Reissner-Nordstrom BH 
solutions. The Lagrangian in this case reads 



(108) 



with fj, > 0. The ESS field solution E{R) is the unique 
positive root of the equation 



2fiE 3 



E 
2~ 



1 



(109) 
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Figure 15. Behaviour of the temperature and specific heat as 
functions of rn and M at constant Q for the RN BH solutions. 
The specific heats exhibit vertical asymptote on the maxima 
of the temperature curves. The slope of the temperature as 
a function of M diverges for EBH configurations. The small 
frames expand the detail of the Cq curves around these EBH 
solutions (see Eqs. (|93[) and (|95[0 . The dashed lines show the 
same functions for the Schwarzschild BHs. 



(R = r/^/Q). As easily seen, this field behaves at the 
center and asymptotically as 



E(R 
E(R - 



0) ~ 2 A1 - 1 / 3 iT 2 / 3 



— > oo 



oo 



R- 



0, 



(110) 



and the model belongs to the A1-B2 family, whose ESS 
solutions in flat space are of finite energy Q. The tem- 
perature of the BH solutions of this model, as a function 
of r/j , can be obtained from Eq. ([68]) in terms of the root 
of ([T09| as 



T(r h ,Q) = J-- 3a /Q- + ^ 



E z 



(111) 



and the expression for the specific heat results directly 
from Eq.(|9"T]>. 

The mass as a function of for several values of Q 
has been plotted in figure |3] for this EH model (main 
frame), where also the behaviour of the RN model is 
displayed (small frame). The only qualitative difference 
between both models corresponds to the small region. 
In the thermodynamically pertinent regions in these fig- 
ures, which lie between the (dotted) lines of extreme 
BHs and the (dashed) line of Schwarzschild BHs, the be- 
haviour of the BH masses versus rh are qualitatively the 
same. Moreover, the behaviours of the functions rj(rh) 
and T(rh) at fixed Q for these models (plotted in figure 
ITTj) as well as the behaviour of the specific heats (plotted 
in figure [15] for the RN model) are also similar. These 
results illustrate the qualitative similarity of the thermo- 
dynamic properties of the BH solutions of UVD and Al 
families near the extreme BH configurations. 
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A first representative example of the A2 family is the 
well known BI model, whose Lagrangian density reads 



1 _ y2 X _ i£y 2 



(112) 



where 1/fi is the maximum field strength. When /i — > 
in (|112l) . the Maxwell Lagrangian (ip = X/2) is recovered. 
The ESS solutions corresponding to model (|112|) take the 
form 



E(R,Q) 



(R = r/^Q). (113) 



The expansion of these solutions around r = 



E(R,Q) = --^ 



3R» 



(R = r/y/Q). (114) 



gives the values of the maximum field strength a = 1/fJ,, 
the coefficients bo = a 3 



co 



= 3a 5 / 128 and the expo- 
nents cr = 4 and A = 8 (see Eo. (j77|) ). The constant A 
in the expansion (|78l) of the Lagrangian density is now 
A = 1/fi 2 = a 2 . 

The expression of the BH temperature as a function of 
Th and Q is obtained from the equations (|68|) . (1112j) and 
(|113[) and reads 



T(r h ,Q) 



1 



47rr^ 
which behaves as 




(115) 



T{r h ^0,Q) 



1 - Q/Qc 



(116) 



for small radius BHs. So we recover the results of the 
precedent analysis on the temperature of the BH config- 
urations close to the extreme {Q = Q c ) and non-extreme 
(Q < Q c ) BPs, described by equation (1541 . For values 
of the charge Q > Q c (A2a cases) the radii of the ex- 
treme BH configurations as functions of Q, obtained from 
Eq.d3U), take the form 



r he (Q > Q c ) = AnQ\ 1 - 



(117) 



and, as expected, the temperature (|115|) vanishes for 
these configurations. 

Similar considerations can be done for the behaviour 
of the specific heat, which results directly from Eqs. ([9T]) . 
(|113|) and (|115[) . For small radius BHs the behaviour of 
Cq is obtained from these formulae or, alternatively, from 
Eq. flMl) in case A2b or from Eq. (|100l) in case A2c, once 



the characteristic parameters of the model have been 
fixed. 

Figure [TBI shows the behaviour of the temperature, for 
these BI black holes, as a function of for several values 
of the charge, corresponding to the cases A2b (two up- 
per curves), A2a (lower curve) and A2c. We have also 
plotted the behaviour of the specific heats for the two 
A2b cases and the A2c case. 
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Figure 16. Behaviour of T(rh,Q) and Cc:(rh,Q) as functions 
of rh for several values of Q corresponding to the different A2 
cases. The picture is obtained from the BI model, as represen- 
tative of the small-rh behaviour of the A2 models with a > 2. 
For the temperature (upper left panel) there are two A2b 
cases (two upper curves, with Q < Q c ), one of them exhibit- 
ing two local extrema (Q < Qc) and the other, obtained for 
lower values of Q, behaving nearly as the Schwarzschild BH 
temperature. The limit = corresponds to the divergent- 
temperature non-extreme BPs. The curve starting at the ori- 
gin corresponds to the critical value of the charge and exhibits 
a EBP there (case A2c). The lower curve, corresponding to 
Q > Qc (A2a cases) exhibits a behaviour similar to that of 
the UVD and Al cases. The associated specific heats are 
plotted in the other panels: upper right and lower left panels 
correspond to the two A2b temperature curves. The lower 
right panel corresponds to the cases A2c, with the curve start- 
ing at ru = 0, but it is also representative of the cases A2a 
if the curve starts on a finite value of ru, corresponding to 
the EBH configurations. The vertical asymptote of Cq cor- 
respond to local extrema of the temperature. 



The A2 models with a = 1 lead to new thermody- 
namic behaviours. In particular, the temperature of the 
extreme black points becomes now finite. Let us give two 
examples of these cases. The first one corresponds to a 
family of NEDs defined through the expression of their 
ESS solutions, which take the form 



E(R,Q) = 



b R + c R > 
1 + i?" 



(118) 



(R = rj^fQ) where the exponents A > 1 and fi > A + 1, 
and the positive constants a, bo and Co are to be cho- 
sen in such a way that E(R) be positive and monoton- 
ically decreasing everywhere, as required for admissibil- 
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ity. From Eq. (|118l) and the first integral ([6]) we can ob- 
tain, through a quadrature, the expression of the admis- 
sible Lagrangian densities supporting these soliton-like 
elementary solutions. We shall make the choice a = 1, 
bo = 0.2, Co = 0.1, \i = 9/2 and two different choices 
for the exponent A (A = 3/2, corresponding to a B3 
asymptotic behaviour or A = 5/2, corresponding to a 
B2 asymptotic behaviour). Both choices lead to admis- 
sible models, for which the coefficients of the expansion 
of E(R) around R = coincide with those of the numer- 
ator of (|118[) . The two choices for the exponent A lead 
to the same value of the extreme BP temperature but to 
different behaviours for the temperature of small-radius 
BHs (see table [TTJ) . For A = 3/2 < 2 the temperature 
decreases from its BP value (= 2bo/y/na) with infinite 
negative slope whereas for A = 5/2 > 2 the temperature 
increases from the same value at rh = with finite slope 
(= 2A). In both cases the specific heats of the extreme 
BP configurations vanish, as 

C Q (r h -+ 0,Q C ) ~ -^Jrf = -2^r][\ (119) 
for the model with A = 3/2, and as 

C Q (r h -+ 0, Q e ) ~ sjl 2 -^r h = ^r h , (120) 

for the model with A = 5/2. 

Figure [17] shows the behaviour of the temperature as 
a function of at fixed Q = Q c for both models. Even 
though this behaviour is similar for finite horizon radius, 
the qualitative difference around = (extreme BP 
configurations) is apparent. 

Finally, let us consider a family of models representa- 
tive of the subclass of the A2 NEDs satisfying the con- 
dition a = 1, A = 3 and A = 0. The assumed form of 
the Lagrangian for X > is 
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Figure 17. Behaviour of T(rh,Q) at constant Q = Q c and 
for small radius BHs associated to the fields (|118[) for two 
choices of the parameter A = 5/2 (upper curve) and A = 3/2 
(lower curve) with fixed values for the other parameters (see 
the text). Note the qualitative difference in the behaviour 
of the EBPs, which exhibit both the same finite temperature 
but with a finite positive slope for A = 5/2 and a negative 
divergent slope for A = 3/2. The small frame expands the 
region of configurations close to the EBPs. 



using the first integral (jSJ) and reduce to solve a quartic 
algebraic equation. Nevertheless, we are interested here 
in the behaviour of the field and the thermodynamic func- 
tions near the extreme BP configurations. By replacing 
the expression (fTTf and the derivative of (|12ip in the first 
integral we determine straightforwardly the parameters 
of the first terms of the polynomial expansion of the field 
around the center as 



a = a ; b = 1 ; c = ; c = 1 ; A = 3, (123) 

or, equivalently, 



<p{X,Y = 0) 



a 



-X 3/2 - X 1/2 



(121) 

where a is a positive constant parameterizing the family. 
It is easy to verify that, as a function of X, at Y = 0, 
this expression increases monotonically from y>(0, 0) = 
and exhibits a vertical asymptote at X = a. Its slope at 
the origin takes the finite value 



d<p 
OX 



x->-o 



(122) 



and, consequently, this is a one-parameter family of ad- 
missible models belonging to the class A2-B2 (obviously, 
for X < these Lagrangians should be properly extended 
taking into account the admissibility conditions). The 
general ESS solution of these models can be obtained ex- 
plicitly in terms of the parameter a and the charge Q by 



E(r h ^0,Q)~a 



(124) 



(R h 

— T^tr) an d we see tnat a must be interpreted as 
the maximum field strength. Moreover, using Eq.([8Tj) 
we are led to the value A = for the integration con- 
stant of the expansion of the Lagrangian around rh = 0. 
Consequently, a glance to tablelTllshows that the temper- 
ature of the extreme BP configuration in these cases is 
Tebp = 2/v / 7ra, the slope of the temperature at 77, = 
vanishes and becomes negative as rh increases, and the 
specific heats attain, for the extreme BP configurations, 
the negative value Cq-ebp = —a 2 /AO. 

Figure [18] displays these temperature and specific heat 
behaviours for the critical (Q — Q c ) configurations, 
which become atypical in the limit of extreme BP con- 
figurations. 
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Figure 18. Behaviour of T(rh,Qc) an d Cq(rh,Qc) for the 
model Q12ip . As in the case of figure fTTl the temperature of 
the EBP configuration is finite, but now its slope vanishes 
there. 



V. FINITE THERMODYNAMIC RELATIONS, 
SCALING GROUP AND GENERATING 
EQUATIONS 

Let us consider now some finite relations arising be- 
tween the state functions introduced in the precedent 
section. 



A. Generalized Smarr law 

In recalculating the derivative of e ex with respect to 
the charge at constant in equation we can use the 
scaling law of the external energy function (|20[) and we 
are led to the relation 



Q& + TS 



3M - r h 



(125) 



or, alternatively, 



M 




\{TS + Q<S>), 



(126) 



where the dependence in has been replaced by depen- 
dence in the entropy S. In the particular case of the RN 
black hole this expression boils down to 



V(r h , Q) = Q<S>(r h , Q) - 2e ex (r h , Q) = 0, (129) 

arises only for the ESS black hole solutions of the 
Einstein-Maxwell model. This condition is related to the 
correspondence between the energy content of the field 
created by a system of charges at rest and their electro- 
static potential e nerg y, which holds in flat-space Maxwell 
electrodynamics [35j, but not in general NEDs. Thus, 
for any admissible G-NED the function V(rh,Q) gives 
the deviations from the exact Smarr law induced by the 
effect of the nonlinear self-couplings of the electric field. 
Owing to the scaling properties (|5Tj) and (|20p this func- 
tion scales as the external energy 



V(r h ,Q) = Q 3 ^V(R h ,i 



(130) 



(Rh = fh/y/Q) aud is fully determined by its character- 
istic expression for Q = 1. 

We conclude that Eq. (1126p is the natural generaliza- 
tion for admissible G-NEDs of the Smarr formula (|127p 
of the Einstein-Maxwell theory. 



B. Other finite relations 

As already mentioned, the thermodynamical approach 
regards the ESS black hole solutions obtained from a 
given gravitating NED as different states of a unique 
system characterized by two state variables, the most 
immediate ones being the integration constants M and 
Q, which arise in the resolution of the Einstein equations 
(1241) . On the other hand, in solving the thermodynamic 
problem for a given model, we start with the first in- 
tegral of the electromagnetic field equations ((6]) and we 
can obtain the state variables as functions of and Q. 
New state variables can be defined as a = a(r^,Q) and 
P — /3(rh,Q), where these functions are assumed to be 
invertible in some pertinent ranges, leading to 



M = 2TS + <2$, (127) 

which is the well known Smarr formula [34j . As another 
example, for the NED family defined in Ref.[18] by the 
Lagrangian density (f(X, Y = 0) = X p (particularized 
here to three space dimensions) , p being an integer num- 
ber, Eq. (|126p provides the result 

M = 2TS +— , (128) 
P 

which is precisely the Smarr-like law obtained in that 
reference. 

Moreover, it is straightforward to show that the equiv- 
alence of (I126P and (|127l) . which requires 



r h = f(a,/3) ; Q = g(a,P), (131) 

where the functions / and g can be multivaluated. 

The universal scaling laws obtained in section ITVl allow 
to build new finite relations between different thermody- 
namic variables, whose form is now dependent on the 
particular model. Let us write these scaling laws in a 
new (but equivalent) form for the main state variables 
introduced so far. First of all, the scaling law (O for the 
electrostatic field leads to 

E{9r hl 6 2 Q) = E(r h ,Q), (132) 

where 9 is any positive parameter (which can be a state 
variable). From Eq. (|75p we obtain for the temperature 
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T(9r h ,9 2 Q) 



1 



4n9r h 

From Eq. (|66|l we obtain for the mass 



9T(r h ,Q). 



(133) 



M(9r h , 9 2 Q) = g(1 2 ° K h + 9 3 M(r h , Q), (134) 

whereas for the electrostatic potential on the horizon we 
obtain, from Eq. ([5"Tj) . 



H9r h ,9 2 Q) =e$(r h ,Q). 



(135) 



Let us denote as T(9) (9 > 0) these scaling transforma- 
tions. From Eqs. (|132j) - ([135|) it is easy to show that the 
product of two transformations satisfies 



L )*r( 



r(0i • 9 2 



(136) 



Consequently the set of scaling transformations for the 
thermodynamical state functions associated to the BH 
solutions of any admissible G-NED has the structure of 
a multiplicative one-parameter group. 

From these laws and the identification of 9 with differ- 
ent state variables we can obtain several finite relations 
whose generic forms are the same for all models, but in- 
volve particular expressions associated to each NED. Let 
us give two simple examples. Consider first the func- 
tion $(rh,Q) obtained from the first integral (0 (after 
a quadrature) for a particular model. Using the scaling 
law (|1 35[) with ff = Twe are led to 



T$ = $ T 



QT 2 



(137) 



where has been written in terms of the entropy. In 
the same way, the explicit expression of the temperature 
(f6"T)) and the scaling law (|1 33[) . with 9 ~ <E>, lead to 



lead to finite relations involving the derivatives of the 
state functions and are also "universal", in the sense 
that they are first-order differential equations whose form 
is independent of the particular gravitating NED. The 
characteristic curves of these equations generate, through 
properly chosen boundary conditions, surfaces in the 
space of state variables involved in the particular scaling 
law. The points of each one of these surfaces correspond 
to BH states associated to a given gravitating NED. As 
the structure of the beam of characteristic curves is also 
universal, this procedure allows for a separate analysis of 
some mo del- independent thermodynamic properties and 
those mo del- dependent behaviours induced by the par- 
ticular boundary conditions corresponding to each NED. 
Let us present here a particular example illustrating the 
method and some of the expected results. 

Using the expression (1521) of the Q-derivative of the 
electrostatic potential at constant Th and replacing the 
electrostatic field on the horizon in terms of the gradient 
of $ we are led to 



= $ 



'V, 



dr h 



(139) 



This equation (which can also be obtained directly by 
derivation of Eq. (|135|) with respect to the parameter 9 at 
6 = 1) involves the function $(r/j, Q) and its first deriva- 
tives with respect to these arguments. It is remarkable 
because its "universality" , since it must be satisfied by 
the solutions of any G-NED, as a direct consequence of 
the scaling law (|135|) . which is a universal property for 
field theories associated to the action (|2"2"|) . On the other 
hand, it is independent of the general laws of thermody- 
namics, as can be easily checked for other thermodynami- 
cal systems. Thus, the particular forms of the "equations 
of state" of the systems of ESS black hole solutions as- 
sociated to the different particular models, characterized 
by the two independent parameters Th and Q, are solu- 
tions of this equation, corresponding to surfaces in the 
$ — r/j — Q space. In order to explore further this inter- 
pretation, let us write Eq. (|139p under the form 



r$ = r $ 



$2 



1 



4$VttS 



(138) 



These equations, together with the generalized Smarr law 
(|126l) . allow for the elimination of two of the five state 
variables &,Q,M,T and S, leading to different expres- 
sions for the "equation of state" of the set of BHs asso- 
ciated to a given model. 



2ga$ 



2S d$ 
IF dS 



-1 = 0, 



(140) 



where the horizon radius has been replaced by the en- 
tropy as independent variable. This expression suggests 
the introduction of the new variables 



x = ln(Q), t/ = ln(S), z = ln($ 2 ) 
in terms of which the equation takes the form 



(141) 



C. Generating equations 

The universality of the scaling laws allows for a more 
direct and enlightening procedure in generating the equa- 
tions of state of general BH configurations. These laws 



dz dz ^ 
dx dy 



(142) 



This is a very simple first-order, linear, partial differen- 
tial equation whose set of characteristics is the beam of 
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straight lines parallel to the vector (1, 1, 1) in the (x, y, z) 
space. Consequently, the equations of state of the sys- 
tem of ESS black hole solutions associated to the G-NEDs 
arc represented in this space by ruled surfaces generated 
by lines of this beam. In order to choose the particu- 
lar surface associated to the system of BH solutions of 
a particular G-NED we can solve, for example, the first 
integral © for the characteristic (Q = 1) case and de- 
termine the form of $(r, Q = 1) through a quadrature. 
This defines a curve z — z(x = 0, y) in the y — z plane, 
which is a boundary condition for Eq. ([142j) (see figure [T5)) 
determining a unique surface solution z = z(x,y) in the 
(x, y, z) space, generated by the beam of characteristics. 
In the ($, Q, S) space the parametric equations for the 
beam of characteristics of (|140|) read 



dM dM 3, r e y' 2 n , 1<lr , 

^ + ^-2 M + 7^ = ' (145) 

whose set of solutions contains the state surfaces M = 
M(x,y) associated with all the admissible G-NEDs. 

These preliminary results require a more detailed in- 
vestigation of the structures which lie at the ground of 
these scaling symmetries and their consequences for the 
BH thermodynamics. A systematic exploration of this 
issue is in progress, but goes beyond the scope of this 
paper. 

VI. CONCLUSIONS AND PERSPECTIVES 




Figure 19. z = ln($ 2 ) versus y = ln(S') curves in Eq. (|141[) 
for the RN solution (dashed line) as compared to the BI and 
EH solutions, in solid style for some choices of the parameters 
and unit charge. 



Q = Q t ; S = S t ; $ = $ V^, 



(143) 



where r is the (positive) parameter and Qo, Sq and $o are 
the coordinates of the points defining the different curves 
of the beam, corresponding to a particular BH configu- 
ration. These curves are parabolas having the origin as 
a common vertex, where the beam is singular. 

A similar analysis can be carried out from the scaling 
laws of other state variables. For example, the scaling 
law for the temperature (|133l) leads to the generating 
equation 



-y/2 



8T dT T 



dx dy 2 y/n 



= 0, 



(144) 



where x and y are the same independent variables de- 
fined in Eq. (|141j ). This is a linear, inhomogeneous, first- 
order, partial differential equation, generating the state 
surfaces T — T(x,y) for the different G-NEDs, once the 
corresponding boundary conditions are specified. In the 
same way, the scaling law for the BH mass (|134l) leads to 
the generating equation 



In this paper we have performed a generalization of 
the well known thermodynamic analysis of the RN black 
hole solutions of the Einstein Maxwell field equations to 
the BH solutions of any admissible gravitating NED. The 
exhaustive classification of the admissible NEDs in terms 
of the central and asymptotic behaviours of their elemen- 
tary solutions in flat space, and the analysis of the geo- 
metric structure of their gravitating counterparts, have 
been carried out in references [l(| and [11(. There, the 
asymptotic flatness of any gravitating ESS solution of 
these admissible NEDs has been established, regardless 
of the asymptotic behaviour of their integral of energy in 
flat space. When this integral diverges at large r (always 
slowly than r, owing to the admissibility conditions) we 
are led to the families of asymptotically IRD models, 
whose gravitating BH solutions exhibit geometric struc- 
tures similar to those of the other models at finite ra- 
dius, but approach flatness at large r slower than the 
Schwarzschild field. For such asymptotically anomalous 
models the standard thermodynamic analysis is meaning- 
less, at least in the usual way (undefined ADM masses 
and asymptotic divergence of the electrostatic poten- 
tials). 

For models whose ESS solutions have asymptotically 
convergent integrals of energy in flat space (B mod- 
els) the gravitating elementary solutions behave as the 
Schwarzschild field at large r, no matter the asymptotic 
degree of convergence of these integrals. In such cases the 
ADM mass of the BH solutions is well defined, the elec- 
trostatic potential can be set to zero asymptotically and 
the thermodynamic laws can be consistently established. 
The thermodynamic properties of the BH solutions in 
these cases have been analyzed in terms of the behaviours 
of the integral of energy and the field around the center 
of the flat space ESS solutions. When this integral (and, 
hence, the field) diverge as r — > (UVD models) the 
geometric structures (two-horizons BHs, single-horizon 
extreme BHs or naked singularities) and the associated 
thermodynamical behaviours are qualitatively similar to 
those of the Reissner Nordstrom BHs, although the tem- 
perature and the specific heats as functions of the hori- 
zon radius (or of the mass) can exhibit more involved 
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behaviours, related to the more or less involved struc- 
ture of E(r) in each case. For models exhibiting finite- 
energy ESS solutions in flat space and whose electrostatic 
field diverges at the center (slower than r -1 , Al mod- 
els) the associated BHs can exhibit new geometric fea- 
tures, as compared with the UVD cases (mainly the ex- 
istence of non-extreme single- horizon BHs), but the form 
of the mass-radius relations and the behaviours outside 
the event horizon are qualitatively similar in both cases 
and hence we are led to similar thermodynamic proper- 
ties. 

For models whose ESS solutions in flat space are finite- 
energy with finite maximum field strength (A2 family, to 
which the Born-Infeld model belongs) the set of associ- 
ated gravitating ESS solutions exhibits a richer structure, 
including extreme and non-extreme black points as new 
possible configurations. These new configurations have 
thermodynamic properties which are related to the be- 
haviour of the electrostatic field around the center. This 
relation has been extensively analyzed and, in particu- 
lar, we have determined the conditions to be satisfied by 
these A2 models in order to support extreme black point 
solutions with zero, finite or divergent temperature. Con- 
cerning the non-extreme black points, their temperature 
is always divergent. 

We have performed a detailed analysis of the proper- 
ties of the extremal BH solutions for the different fam- 
ilies. The thermodynamic behaviour of these configu- 
rations differ strongly among the different families and 
suggest the possible existence of models exhibiting van 
der Waals-like first-order phase transitions in the $ — Q 
plane. 

Once the possible structures and thermodynamic prop- 
erties of the BH solutions of the different families of ad- 
missible NEDs have been determined and classified, the 
results of the present analysis allow for the explicit con- 
struction of specific models whose solutions exhibit pre- 
scribed behaviours. In this way we have built some new 
models as examples illustrating different aspects of the 
study carried out here. 

From the invariance of scale of the general NED field 
equations in flat space, we have established universal scal- 
ing laws, satisfied by the thermodynamic functions defin- 
ing the BH states of any model. These scaling laws lead 
to universal first-order partial differential equations in- 
volving sets of three state variables, whose solutions are 
surfaces in the three-dimensional spaces defined by these 



variables. Each one of these surfaces is generated by 
the beam of characteristics associated to the equation. 
From appropriate boundary conditions we can select the 
surface solution associated to a particular G-NED. The 
points of this surface correspond to the BH states which 
are the solutions of this model, and its analytic equation, 
relating the involved variables, defines a expression of the 
"equation of state" of this system of BHs. 

To conclude, let us mention some perspectives and pos- 
sible extensions of this study: 

First, in Rcf.|8] we have established that the analysis 
of the ESS problem for NEDs in flat space can be ex- 
tended to the same problem for generalized nonabelian 
gauge field theories of compact semi-simple Lie groups. 
More explicitly, we have shown that the ESS solutions of 
a generalized gauge field theory, defined by a Lagrangian 
density of the form (fTJ) (in terms of the usual gauge invari- 
ants with the ordinary trace definition in this nonabelian 
case) have the same functional dependence in r as those 
of the NED defined by a Lagrangian density of the same 
form (in terms of the abelian gauge invariants) . In this 
way we have recently communicated some preliminary 
results on the extension of the analysis performed here 
to systems of nonrotating BHs of nonabelian origin [36| . 
A more detailed report on this issue is in progress. 

Second, the extensions of the present general analysis 
(geometric structures and thermodynamics) to the cases 
of NEDs coupled to more general actions for the grav- 
itational field (e.g. with a cosmological constant term, 
or Gauss-Bonnet and Lovelock theories) is in progress, a 
context where some particular cases of NEDs have been 
already explored in the literature [2r| . 

Finally, let us stress that a more complete analysis of 
the scaling symmetries and their consequences on the BH 
thermodynamics is necessary. In particular, the study 
of the structures of the beams of characteristics associ- 
ated to different generating equations and of the bound- 
ary conditions selecting the equations of state of the BH 
systems associated with different NEDs, allows for the 
exploration of eventual phase transitions, an issue also 
under development. 
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